Accepted for The Astrophysical Journal 



Fermi observations of TeV-selected AGN 

A. A. Abdo^'^ M. Ackermann^ M. Ajello'', W. B. Atwood^ M. Axelsson'^'^ L. Baldini^ 
J. Ballet^ G. BarbielliniiO'", M. G. Baring^^, D. Bastieri^^-i^^ K. Bechtol^ R. Bellazzini^ 

B. Berenji^, E. D. Bloom^, E. Bonamente^^'^^, A. W. Borgland^, J. Bregeon^, A. Brez^, 
M. Brigidal^'l^ P. Bruel^^ T. H. Burnett^o, G. A. Caliandrol^'l^ R. A. Cameron^ 
R A. Caraveo2\ J. M. Casandjian^ E. CavazzutP, C. Cecchi^^.ie^ q ggiji^23,24,25^ 

A. Chekhtman2'26, C. C. Cheung^^, J. Chiang^ S. CipriniiS'ie^ R. Claus^ 
J. Cohen-Tanugi^^, L. R. Cominsky^^, J. Conrad2^'^'^°'^\ S. Cutini^^, A. de Angelis^^, 
F. de Palma^^•l^ G. Di Bernardo®, E. do Couto e Silva^ R S. Drell^ A. Drlica-Wagner^ 
R. Dubois'^, D. Dumora^^'^'', C. Farnier^"^, C. Favuzzi^'''^®, S. J. Fegan^^'\ J. Finke^'^, 
W. B. Focke^ P. Fortin^^ L. Foschini^s, M. Frailis^^, Y. Fukazawa3^ S. Funk^ 
P. Fusco^'^'^^ F. Gargano^®, D. Gasparrini^^ N. Gehrels^^'^"^, S. Germani^^'^^ G. Giavitto^^ 
B. Giebels^^ N. Gigliettol^•l^ P. Giommi^^, F. Giordanol^'l^ T. Glanzman^ G. Godfrey^ 

I. A. Grenier^, M.-H. Grondin^^'^"^, J. E. Grove^, L. Guillemot^^'^'^, S. Guiriec^^, 
Y. Hanabata3^ M. Hayashida^, E. Hays^^ D. Horan^^ R. E. Hughes^", M. S. Jackson^^, 

G. Johannesson^ , A. S. Johnson^, R. P. Johnson^, W. N. Johnson^, T. Kamae^, 
H. Katagiri^e, J. Kataoka^^'^^^ N. Kawai^l'^^ M. Kerr^o, J. Kn5dlseder^^ M. L. Kocian^ 
M. Kuss®, J. Lande^, L. Latronico®, M. Lemoine-Goumard^^'^'^, F. Longo^°'^^, 
F. Loparcol^•l^ B. Lott^^'^^, M. N. Lovellette^, P. Lubranol5'l^ G. M. Madejski^ 
A. Makeev2'26, M. N. Mazziotta^^ W. McConville^^'S^^ J. E. McEnery^^, C. Meurer^^'^, 
P. F. Michelson^ W. Mitthumsiri"^, T. Mizuno^^ A. A. Moiseev^^'^"^, C. Monte^^'^^ 
M. E. Monzani^ A. Morselli^^ 1. V. Moskalenko^ S. Murgia^ P. L. Nolan^ J. P. Norris^^ 

E. Nuss^^ T. OhsugP, N. Omodei^ E. Orlando^^ J. F. Ormes^^ M. Ozaki^^ 
D. Paneque^, J. H. Panetta^, D. Parent^^'^^, V. Pelassa^^ M. Pepe^^'^^ M. Pesce-Rollins^ 
F. Piron^^, T. A. Porter^ S. Raino^^•l^ R. Rando^^'^^, M. Razzano^ A. Reimer^^-^, 

O. Reimer^9,4^ t. Reposeur33'34^ l C. Reyes^°, S. Ritz^ L. S. Rochester^ 
A. Y. Rodriguez^i, M. Roth^o, F. Ryde^O'^ H. F.-W. Sadrozinski^ D. Sanchez^^.i^ 
A. Sander^°, P. M. Saz Parkinson^ J. D. Scargle^^ T. L. Schalk^ A. Seller hoW^'^, 
C. Sgr6^ M. S. Shaw^ E. J. Siskind^^ D. A. Smitli^^-s^^ P. D. Smith^°, G. Spandre^ 
P. Spinellil^'l^ M. S. Strickman^, D. J. Suson5^ H. Tajima^ H. TakahashP, 
T. Takaliashi^® , T. Tanaka^, Y. Tanaka^®, J. B. Thayer^, J. G. Thayer^, D. J. Tliompson^^, 
L. Tibaldol3'9'^^ D. F. Torres55'5\ G. Tosti^^'^^ A. Tramacere^-^^ Y. Uclliyama^S'^ 
T. L. Usher^ V. Vasileiou23'24,25^ n. Vilchez^^ V. Vitale^^-^^ A. P. Waite^ P. Wang^ 
B. L. Winer^°, K. S. Wood^, T. Ylinen^o-ss,?^ Ziegler^ 



- 2 - 



^Corresponding authors: S. J. Fegan, sfegan@llr.in2p3.fr; D. Sanchez, dsanchez@llr.in2p3.fr. 

^Space Science Division, Naval Research Laboratory, Washington, DC 20375, USA 

^National Research Council Research Associate, National Academy of Sciences, Washington, DC 20001, 
USA 

*W. W. Hansen Experimental Physics Laboratory, Kavli Institute for Particle Astrophysics and Cosmol- 
ogy, Department of Physics and SLAC National Accelerator Laboratory, Stanford University, Stanford, CA 
94305, USA 

^Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and 
Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064, USA 

^Department of Astronomy, Stockholm University, SE-106 91 Stockholm, Sweden 

'^The Oskar Klein Centre for Cosmo Particle Physics, AlbaNova, SE-106 91 Stockholm, Sweden 

^Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, 1-56127 Pisa, Italy 

"Laboratoire AIM, CEA-IRFU/CNRS/Universite Paris Diderot, Service d'Astrophysique, CEA Saclay, 
91191 Gif sur Yvette, France 

^°Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 1-34127 Trieste, Italy 

^^Dipartimento di Fisica, Universita di Trieste, 1-34127 Trieste, Italy 

i^Rice University, Department of Physics and Astronomy, MS-108, P. O. Box 1892, Houston, TX 77251, 

USA 

^•^Istituto Nazionale di Fisica Nucleare, Sezione di Padova, 1-35131 Padova, Italy 
^■^Dipartimento di Fisica "G. Galilei", Universita di Padova, 1-35131 Padova, Italy 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Perugia, 1-06123 Perugia, Italy 
^^Dipartimento di Fisica, Universita degli Studi di Perugia, 1-06123 Perugia, Italy 
^'^Dipartimento di Fisica "M. Merlin" dell'Universita e del Politecnico di Bari, 1-70126 Bari, Italy 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Bari, 70126 Bari, Italy 
^^Laboratoire Leprince-Ringuet, Ecole polytechnique, CNRS/IN2P3, Palaiseau, France 
2°Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 
^^INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica, 1-20133 Milano, Italy 
^^Agenzia Spaziale Italiana (ASI) Science Data Center, 1-00044 Frascati (Roma), Italy 
23NASA Goddard Space Flight Center, Grecnbelt, MD 20771, USA 

^''Center for Research and Exploration in Space Science and Technology (CRESST), NASA Goddard Space 
Flight Center, Greenbelt, MD 20771, USA 

^^University of Maryland, Baltimore County, Baltimore, MD 21250, USA 



-3- 



^^George Mason University, Fairfax, VA 22030, USA 

^''Laboratoire de Physique Theorique et Astroparticules, Universite Montpellier 2, CNRS/IN2P3, Mont- 
pellier. Prance 

^^Department of Physics and Astronomy, Sonoma State University, Rohnert Park, CA 94928-3609, USA 

^^Department of Physics, Stockholm University, AlbaNova, SE-106 91 Stockholm, Sweden 

3°Department of Physics, Royal Institute of Technology (KTH), AlbaNova, SE-106 91 Stockholm, Sweden 

^^Royal Swedish Academy of Sciences Research Fellow, funded by a grant from the K. A. Wallenberg 
Foundation 

^^Dipartimento di Fisica, Universita di Udine and Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, 
Gruppo CoUegato di Udine, 1-33100 Udine, Italy 

^•^Universite de Bordeaux, Centre d'Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33175, 
France 

3^CNRS/IN2P3, Centre d'Etudes Nucleaires Bordeaux Gradignan, UMR 5797, Gradignan, 33175, France 
^^INAF Osservatorio Astronomico di Brera, 1-23807 Merate, Italy 

^^Department of Physical Sciences, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 
^^University of Maryland, College Park, MD 20742, USA 

^^Istituto Nazionale di Fisica Nucleare, Sezione di Trieste, and Universita di Trieste, 1-34127 Trieste, Italy 
^^University of Alabama in Huntsville, Huntsville, AL 35899, USA 

^•^Department of Physics, Center for Cosmology and Astro-Particle Physics, The Ohio State University, 
Columbus, OH 43210, USA 

^^Department of Physics, Tokyo Institute of Technology, Meguro City, Tokyo 152-8551, Japan 

^^Waseda University, 1-104 Totsukamachi, Shinjuku-ku, Tokyo, 169-8050, Japan 

^^Cosmic Radiation Laboratory, Institute of Physical and Chemical Research (RIKEN), Wako, Saitama 
351-0198, Japan 

^^Centre d'Etude Spatiale des Rayonnements, CNRS/UPS, BP 44346, F-30128 Toulouse Cedex 4, France 
^^Istituto Nazionale di Fisica Nucleare, Sezione di Roma "Tor Vergata" , 1-00133 Roma, Italy 

Department of Physics and Astronomy, University of Denver, Denver, CO 80208, USA 
^''Max-Planck Institut fiir extraterrestrische Physik, 85748 Garching, Germany 

^^Institute of Space and Astronautical Science, JAXA, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, 
Japan 

^^Institut fiir Astro- und Teilchenphysik and Institut fiir Theoretische Physik, Leopold-Franzens- 
Universitat Innsbruck, A-6020 Innsbruck, Austria 

^°Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA 



-4- 



ABSTRACT 



We report on observations of TeV-selected AGN made during the first 5.5 
months of observations with the Large Area Telescope (LAT) on-board the Fermi 
Gamma-ray Space Telescope (Fermi). In total, 96 AGN were selected for study, 
each being either (i) a source detected at TeV energies (28 sources) or (ii) an 
object that has been studied with TeV instruments and for which an upper-limit 
has been reported (68 objects). The Fermi observations show clear detections of 
38 of these TeV-selected objects, of which 21 are joint GeV-TeV sources and 29 
were not in the third EGRET catalog. For each of the 38 Fermi -detected sources, 
spectra and light curves are presented. Most can be described with a power law 
of spectral index harder than 2.0, with a spectral break generally required to 
accommodate the TeV measurements. Based on an extrapolation of the Fermi 
spectrum, we identify sources, not previously detected at TeV energies, which 
are promising targets for TeV instruments. Evidence for systematic evolution 
of the 7-ray spectrum with redshift is presented and discussed in the context of 
interaction with the EBL. 

Subject headings: Gamma rays: observations 



1. Introduction 



At energies above approximately 100 GeV (hereafter the TeV energy regime), ground- 
based 7-ray observatories have detected 96 sources over the past two decades. The pace of 
discovery in this energy regime has been particularly high since the inceptio n of the latest 
generation of instruments: H.E.S.S., CANGAROO, MAGIC and VERITAS dWeekesI I2OO8I . 
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for recent review). Online catalogs, such as TeVCafl present continuously updated views of 
the TeV 7-ray sky. The majority of the TeV sources are galactic, however 30 extragalactic 
sources have also been detected, of which 28 correspond to Active Galactic Nuclei (AGN), the 
other two being recently detected starburst galaxies. Most (25) of these TeV AGN are blazars, 
an AGN sub-category in which the jet of relativistic plasma ejected from the core is roughly 
coaligned with our line-of-sight and hence appear Doppler boosted. The majority (24) of 
the TeV blazars belong to a further sub-category, the BL Lac objects (from BL Lacertae, 
the prototype for the class), which do not have significant emission or absorption features in 
their optical spectra, making it difficult to measure their redshift directly. 



The first blazar detected at TeV energies was Markarian 421 (Mrk 421; iPunch et al. 



19921 ). at a redshift of z = 0.031. It is seen to be highly variable, with fiux varying between 
~ 0.15 and > 10 the fiux of the Crab N ebulaE ((t>crah)- Dou bling timescales as short as 15 
minutes have been observed during fiares (IGaidos et al.lll996l ). Mrk 421 has a hard spectrum, 
with mean photon inde x of F = 2.5, and has shown clear evidence of spectral variability 
during fiaring episodes (IKrennrich et al.ll2002al ). The detection of such a distant object was 
interp reted in the context of 7-ray attenuation through pair-production (iGould fc Schreder 
19671) to produce a limi t on the power density of extragalactic background light (EBL; 
Stecker fc de Jagerl 119931). Mrk 421 was detected by EGRET, the predecessor of the Fermi 



Large Area Telescope (IThompson et al 



third EGRET catalog (3EG; iHartman et al. 



1993 



Kanbach et al.lll9881) . and was reported in the 
1999h with a detection significance of ~ lOa. 



However, EGRET did not have sufficient sensitivity to make detailed measurements on the 
short timescales required to match the TeV observations. 

Since this initial discovery, 28 AGN sources have been detected at TeV energies, the 
most distant reported being 3C 279, at a redshift of z = 0.54. Like Mrk 421, many of 
these have shown evidence for variability, undergoing episodic fiaring activity with short 
doubling timescales. To date, the rnost ex treme example of variability has been observed 
from PKS 2155-304 (lAharonian et al.ll2007al ). in which the fiux was seen to reach a maximum 
of ~ 15 0crab with doubling times as short as 225 s. However, a pproximately half of the TeV 
blazars show n o evidence for variabili t y, e.g. PKS 2005-4 89 (lAharonian et al.l l2005al ) and 
PG 1553-1-113 (lAharonian et al.ll2006d : lAlbert et al.ll2009l ). The detections of more distant 
objects lead to tighter constraints on the the level of the E BL, suggesting that its d ensity 
is close to the minimum level required from galaxy counts (lAharonian et al.ll2006al ). The 



^http://tevcat. uchicago.edu, see Wakelv fc HoranI ( 2008 ) 

^The Crab Nebula is the brightest steady TeV source and the "standard candle" of TeV astronomy, 
defining the "Crab Unit" in which TeV fluxes and limits ar e often expressed. We adopt the value of ^crab(> 



E) = 2.1 X 10-"(£:/lTeV)-i-5cm 



from 



Hillas et al 
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spectra of the majority of TeV blazars are adequately described by a simple power with 
r > 2, and with those of the more distant source s being considerably softer, up to F ^ 4 (e.g. 
Acciari et al1l2009bl : lAlbert et al.ir2007bl . l2008bl : kharonian et al.l 120063 . boOSah . The peak 
in the measured 7-ray spectrum (in uF^, representation) of these objects lies outside of the 
energy range of the TeV instruments, making it difficult to fully constrain models of emission 
using TeV observations alone. The spatial and spectral properties of the TeV-detected AGN 
are presented in Tables [1] and [21 together with references to TeV observations of each source. 

The Large Area Telescope (LAT) is a pair conversion telescope on the Fermi Gamma- 
ray Space Telescope ( formerly GLA S T), la unched in June 2008. The Fermi LAT instrument, 
described in detail in lAtwood et al.l (120091 ) . detects 7 rays with energies between 20 MeV and 
>300GeV (hereafter denoted the GeV energy regime). The bulk of the Fermi observational 
program is dedicated to a sky-survey, in which the full 7-ray sky is observed every 3 hours. 
This survey is optimized to produce a uniform exposure to the sky on timescales of months 
and to facilitate the monitoring and detection of variable and flaring 7-ray sources on shorter 
timescales. In its first three months of operation the Fermi LAT mapped the 7-ray sky with 
a sensitivity and precision that exceeds any previous space mission in this energy regime. 
A list of the 205 brightest sources (lOcx or greater) found during that period, and their 
properties, has been published by the Ferrni-LAT collaboration to guide multiwavelength 
studies with other instruments (lAbdo et al.ll2009bl ). This collection of sources is henceforth 
referred to as the OFGL list, with individual sources denoted as OFGL JHHMM.MiDDMM. 
In addition, an in-dept h study of the popu lation of OFGL sources most likely associated with 
AGN has been made (lAbdo et al.ll2009a] ): this population is commonly referred to as the 
LAT bright AGN sample (LBAS). 

Only eight of the 28 TeV-detected AGN were detected by EGRET and included in 
the 3EG catalog (see Table H]); six are BL Lacs. The majority of these 3EG GeV-TeV 
sources were discovered as TeV emitters only with the advent of the latest generation of 
TeV observatories. With the previous generation of instruments only two such extragalactic 
GeV-TeV sources were established, despite dedicated programs to observe 3EG sources (e.g. 
Fegan et al.ll2005l ). indicating the degree of mismatch between the sensitivities and effective 
energy ranges of th ese instruments. In accordance with the blazar sequence hypothesis 
( iFossati et al.lll999l ). BL Lacs are the least luminous class of blazars in the GeV regime, 
but their emission peaks at a higher energy and they have hard photon indices in the GeV 
domain. In fact, above a few 10s of GeV, BL Lacs are often relatively brighter than other 



•^The exceptions being Mrk 421 and Mrk 501 which show evidence for spectral curvature and PKS 2155- 
304 whose spectrum was found to have the form of a broken power-law, during the flaring episode of June 
2006. 
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AGN. With a rapidly falling sensitivity above ~ 5 — 10 GeV, EGRET preferentially detected 
the more luminous, lower-energy-peaked blazars. In contrast, the Fermi LAT has a relatively 
flat effective area at high energies (8000 cm^ for on-axis 7 rays at E > 1 GeV) and an energy 
response that extends beyond 300 GeV, making it much more sensitive to the hard BL Lac 
sources than was EGRET. Of the OFGL sources, 14 are AGN detected at TeV energies. 

One of the most powerful tools for probing the physics underlying the emission from 
AGN is the dedicated multiwavelength observational campaign, in which simultaneous obser- 
vations are made across the full spectrum. Generally, such observations of the TeV blazars 
reveal two non-thermal components: one at lower energies, peaking in the UV to X-ray 
regime, and showing clear evidence of polarization, and a second peaking in the 7-ray regime. 
The low-energy component is commonly interpreted as resulting from synchrotron emission 
from relativistic electrons in the jet, while the high-energy component results from a different 
process, such as inverse-Compton scattering of lower-energy photons in the region of the jet, 
or the decay of 7r° particles produced in interactions of relativistic protons. In many such 
campaigns significant correlation between the X-ray an d TeV 7-ray emission has been de- 



tected (e.g. iBuckley et al.lll996l : lAharonian et al.ll2009bl ). suggesting that a single population 



of rel ativistic particles is responsible for the emission in both regimes (see iKatarzyhski et al. 
2OO5I . for further discussion of X-ray/TeV correlations). Leptonic mechanisms, such as the 



synchrotron self-Compton (SSC) and external- Compton (EC) processes, and hadronic mech- 
anisms have been invoked to explain the broad-band emission and correlated time variability 
seen between the low-energy and high-energy components (see lBottcherll2007l . for further dis- 
cussion). Instances of is olated 7-ray variability have also been detected in some cases (e.g. 
Krawczynski et al.ll2004j ). and it is probable that no simple mechanism will fully explain the 
considerable variety of the blazar behaviour. 

Despite the participation of instruments across a wide range of the spectrum, until 
recently, multiwavelength campaigns have been largely unable to probe the energy range 
between ~ 150 keV and ~ 150 GeV, as no instrument with sensitivity matched to the day- 
to-week timescales of typical campaigns has existed. As such, although the synchrotron 
component has been well measured from radio to X-ray, the full extent of the higher energy 7- 
ray component has not. In the case of the TeV blazars, ground-based 7-ray instruments have 
generally measured the falling edge of the high-energy component (in vFy representation), 
which has usually been consistent with a featureless power-law. The rising edge and peak of 
the emission have, however, been inaccessible, and hence models of the emission have been 
unconstrained in this energy range. In many cases, very different emissions mechanisms 
have been invoked, and can explain the data equally well. With the launch of Fermi, which 
has the sensitivity to measure the emission from the brighter TeV blazars on the day to 
week timescales, a large part of this gap in coverage has been closed. A recent multi- 
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wavelength campaign on PKS 2155-304 was the first in which the r ising and faUing edges o f 



the high-energy SED were simultaneously measured with precision (lAharonian et al.ll2009d ). 
Measurement of the full broadband spectrum and the pattern of correlation between the 
optical, X-ray, GeV 7-ray and TeV 7-ray emission removes degeneracies in modeling of this 
object which were present in the results from previous campaigns that could not measure 
the high-energy component fully. 

In this paper we present the results of the first 5.5 months of Fermi-LAT observations 
of the known TeV blazars and of those AGN for which upper-limits exist at TeV energies. 
The motivation for this study is two-fold: (i) to present as complete a picture of the high- 
energy emission as possible by combining the GeV and TeV results on these objects, and 
(ii) to help guide future TeV observations. For each object detected by Fermi a power-law 
fit to the GeV spectrum is presented, as are light curves on monthly timescales. For the 
brighter sources, lightcurves on 10-day timescale are also given. The GeV power-law spectra 
are extrapolated to TeV energies assuming absorption on the EBL, yielding predictions for 
TeV emission for the simplest case where there is no curvature in the intrinsic spectra of the 
objects. For those objects which are not detected by the Fermi LAT, upper limits in the 
GeV range are presented and extrapolated to TeV energies. 



2. Sources 

The primary objects selected for this study are the 25 blazars and 3 radio galaxies 
detected at TeV energies. These are listed in Table [1] with their coordinates, the AGN 
sub-class of the object, redshift and references to the initial detection at TeV energies. 
In summary, 19 high-frequency-peaked BL Lacs (HBL), 3 intermediate-frequency-peaked 
BL Lacs (IBL), 2 low-frequency-peaked BL Lacs (LBL), 1 flat-spectrum radio quasar (FSRQ) 
and 3 Fanaroff- Riley radio galaxies (type FRl) have been detected by TeV instruments, the 
most distant having a redshift of z = 0.54. Table [2] lists the parameters of a power-law fit 
to the TeV spectra for these objects, where available: the integral flux (0 ± A0) and photon 
index of the fit (F ± AF) and the threshold energy for the observation (i^thres), such that the 
differential spectrum is: 

— = (r-i)^(^y^ = F ( ^ V 

dE ^ ^^thres V^thres; '°°V200GeVy' 

The differential flux at 200 GeV (the median threshold of the measurements), -F200; is cal- 
culated from the TeV power-law spectrum and presented in the table to compare the TeV 
objects at a single energy lying within the domain of the Fermz-LAT observations. For 
some objects, multiple TeV spectra have been measured, either by different instruments, in 
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different epochs, or when the object is in different ffux states. Where possible, the spectrum 
corresponding to a low-ffux state is hsted. 

In addition we search for GeV emission from 68 objects f or which TeV upper-limits were 



published from observatio ns with the Whipple 10m telescope (IHoran et a 



and H.E.S.S. (lAharonian et al.ll2005c 



20031: iFalcone et al.ll2004[). HEGRA (Aharonian et al.l l2004[). MAGIC flAlbert et al. 



12004 



de la Calle Perez et al. 
2008ah 



2008bl ). These targets are listed in Table [3] with the 



lowest flux upper-limit published. 

From these 96 target objects, the 18 lis ted in Table 



sources in the OFGL list and LBAS sample (lAbdo et al 



H are identified or associated with 



2009 



aa. 

Those lists were limited 



to sources with TS > 100 in three months of Fermi-LAT datal^; in this study the criterion to 
claim a detection and derive a spectrum is lowered to TS > 25 and the period of observation 
is increased to 5.5 months. A total of 38 sources are detected by Fermi, of which 21 are 
jointly detected at GeV and TeV energies. We also give an upper-limit for TeV sources not 
detected by Fermi. 



3. Analysis 

Fermi-LAT data from the 5.5 month interval from MJD 54682 to MJD 54842 are pro- 
cessed with the standard analysis chain ScienceTools (ST; version V9R11). The latest 
instrumental response functions (IRFs; version P6_V3) are used to characterize the PSF 
and effective area during the analysis. These IRFs offer a distinct improvement over those 
used in the OFGL analysis, which did not properly account for the presence of remnants of 
non-triggering events in the tracker. This change results in a systematic increase of ~ 15% 
in the derived flux from 7-ray sources, and a possible change in the spectral index, which is 
most pronounced for softer sources. However, despite these improvements, the P6_V3 IRFs 
are based on pre-fiight calibrations, and to be conservative, only eve nts in the eiiergy ra nge 



from 200 MeV to 300 GeV are retained for analysis (see for example. lAbdo et al.ll2009bl ). 



The data for each of the AGN targets is analyzed in a identical manner. Low-level 
processing of the spacecraft data is applied automatically in a pipeline, reconstructing the 
energy, arrival direction and particle type of the primary. Events reconstructed from a 
region of interest (ROI) of 10° around the target location are extracted from this database 
and filtered such that only those having the highest probability of being a photon (those in 
the so-called "diffuse" event class) and having an angle of < 105° with respect to the local 



''The test statistic, TS, is roughly indicative of the significance of the LAT detection of the source squared. 
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zenith (to suppress the background from the Earth albedo) are retained. 



For each target, a background model is const ructed, consisting o f a diffuse galactic 
component, predicted by the GALPROP program Jstrong et al.l l2004al lb[) . a diffuse power- 
law component (for the extragalactic and i nstrumental backg round) and any of the point 
sources from the Fermi three-month catalog ( Abdo et al. 2009b 1^ which overlap the ROI. The 
spectrum for each of the point sources is modeled as a power-law. An unbinned maximum 
likelihood method (ICashlll979l : iMattox et al.lll996l ). implemented as part of the ST by the 
gtlike program, is used to optimize the parameters to best match the observations. 

The validity of the optimized emission model is verified by producing a TS map for 
each region. This is done using the ST gttsmap program, which adds a test source to each 
location over a prescribed region and calculates the improvement in log-Likelihood with 
the inclusion of the test source. Statistically compelling sources, not accounted for by the 
model, are identified visually in the map and added to the background in another iteration 
of gtlike. In addition, high resolution TS maps are produced for each source of interest 
and the centroid of the emission and contours defining the 68%, 95% and 99% probability 
regions calculated. During the construction of the OFGL list, a systematic error of ~ 1' in 
the reconstruction of the centroids of emission from well known bright sources was identified, 
and this error is folded into the emission contours displayed in Figure [H If the centroid of 
emission on the map is Tc and the distance from the centroid to the (statistical) error contours 
are defined parametrically by the functions rl^^^{6), where i = 1,2,3 for the 68%, 95% and 
99% probabilities (P* = 0.68,0.95,0.99), then the contours which account for systematic 
errors are defined as, 



syst 



{0) 



\/rAatW' + (v^21n(l-P0xlT. 



For well detected sources the error contours are roughly circular, with the systematic error 
dominating. For weakly detected sources the contours can be irregularly shaped, with the 
statistical component dominating. 

Among the outputs from the gtlike program are the optimized values of the model 
parameters, the covariance matrix describing their variances and correlations, and the TS of 
each source, which indicates the significance of the source detection. An error contour (called 
butterfly diagram) is computed from these values and plotted to indicate the la confidence 
range of the fitted model. If the power-law is written as F{E) = dN/dE = Fq{E/Eq)~^, 
with normalization parameter Fq ± AFq, photon index F ± AF and covariance cov{FQ,r), 



^We use an internal version of catalog which is not limited to TS>100. 
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the the contour is defined by: 



The narrowest point in the butterfiy occurs at = -Eq exp[cov(Fo, r)/Fo AF^], the so-called 
decorrelation energy. For each source, the butterfiy is drawn between the lowest energy used 
in the analysis (either 0.2 GeV or 1 GeV) and the energy of the highest photon detected from 
the source, subject to the constraint of < 300 GeV. 

For the sources with a detection significance of TS < 25 (~ 5cr), upper limits on the 
integral fiux above 200 MeV are computed assuming a photon index arbitrary fixed at 1.5 and 
2.0. For the very bright LAT sources, the energy range is divided into two bins (200 Me V- 
1 GeV and 1 GeV-300 GeV), and spectra are fitted to each bin separately. This analysis is 
limited to sources for which TS > 100 in each of the energy bins, ensuring that a sufficiently 
accurate spectrum can be derived in each. For all detected sources, systematic errors on 
the fiux and index of the power-law fit to the full Fermi energy range, caused by systematic 
errors in the IRFs used in the analysis, are evaluated using the "bracketing" method of 



Abdo et al.1 (|2009d|). 



In order to make predictions for the TeV energy domain and to make comparisons be- 
tween the Fermi-LAT and TeV spectra, the best-fit spectrum and butterfiy are extrapolated 
up to 10 TeV. This assumes that the intrinsic spectrum of the emission is described by 
a single power-law extending over the GeV and TeV energy range, the simplest and least 
model-dependent assumption that can be made. Above a few hundred GeV, the photons 
interact with the infrared photons from the EBL as they pr opagate through the Uni verse, 



modifying the detected spectrum from a simple power-law. iFranceschini et al.l (120081 ) pro- 
vide tabulated values of the optical depth as a function of the redshift, which are used to 
compute the fiux detectable between 200 GeV and 10 TeV from the extrapolated GeV spec- 
trum. Their EBL model is consistent with experimental measurements: the lower limits 
from galaxy counts and upper limits from observations of TeV blaz ars, and is widely used . 



However it is not necessarily the final word on EBL density (see e.g. iKrennrich et al.ll2008l ). 
and any errors in the model will propagate into the extrapolation of the GeV spectra to TeV 
energies. In addition to absorption o n the EBL, TeV photons may undergo absorption in the 



neighborhood of the source (see e.g. lAharonian et al.ll2008d ). which must also be modeled 



and accounted for to unfold the intrinsic accelerated spectrum from the detected spectrum. 
However, such modehng is beyond the scope of this paper. 

Light curves for each source are produced with time bins of 10 and 28 days (following 
the lunar cycle). The light curves are produced by binning the events by their arrival times 
and performing an independent likelihood analysis for each temporal bin with the same 
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model (same background sources and number of free parameters) as in the fitting of the 
time-averaged spectrum. The probabihty that the hght curve is consistent with being flat, 
from a fit to a constant value, is also computed, and used to evaluate the hypothesis that 
the fitted spectra, averaged over the full period, are valid. 

4. Results and discussion of individual sources 

Of the 28 TeV-selected sources studied here, 21 are detected by Fermi with TS > 25. 
This degree of connection between the TeV blazars and the GeV regime was not found by 
EGRET and the previous generation of TeV instruments, and is evident now only as a result 
of the improved sensitivity and greater overlap between the effective energy ranges of Fermi 
and the current generation of TeV instruments. 

The majority of the TeV blazars detected by Fermi have a photon index F < 2 in the 
GeV regime, the median index is F = 1.9. In contrast, the populations of 42 BL Lacs and 
57 FSRQs from the LBAS sample have median indices of F = 2.0 and F = 2.4 respectively. 
The TeV blazars are amongst the hardest extragalactic objects detected by Fermi. 

Of the 68 extragalactic objects with TeV limits which were considered here, a total of 17 
are detected in the GeV regime. These too have a hard median index, of F = 1.95, indicating 
that they, perhaps, are good targets for deeper observation with TeV instruments. Of these 
17, only one is not a blazar (NGC 1275, an FRl radiogalaxy) , one is an FSRQ (3C 273) and 
the remainder are BL Lacs. That the majority of these 16 blazars are BL Lacs, rather than 
following the ratio of FSRQs to BL Lacs found in LBAS, is most likely a result of the way 
these objects were selected for observation originally by the various TeV groups, rather than 
anything inherently fundamental. These 15 objects break down as follows: 3 LBL, 1 IBL 
and 11 HBL. 

For each of the 38 sources detected by Fermi, power-law fits to the data and extrap- 
olations into the TeV regime are presented in Figure [2] for GeV-TeV sources, and in the 
online material for the GeV sources with TeV upper-limits. Additionally, for TeV sources 
not detected by Fermi, upper limits on the spectra are presented in Figure [31 To justify 
the vahdity of these "averaged spectra" , light curves on 28-day timescales are presented in 
Figure m for all GeV-detected objects, and on 10-day timescales in Figure O for the brighter 
GeV emitters. In the 5.5 months of data analyzed here, the majority of sources do not show 
evidence of flux variation on the timescales tested. For a subset of the sources, as discussed 
below, a TS map is presented in Figure [H The remainder are available in the online material. 

In addition to the figures, the parameters of the power-law spectra and variability are 
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given in tables. For each of the sources detected by Fermi, Table [5] lists 

• the flux and index of the power law, the statistical and systematic errors on these 
quantities, the decorrelation energy, 

• the energy of the highest and fifth highest photons detected within 0.25° of the source 
position, which corresponds to > 99.9% containment according to the P6_V3 IRFs, 
and 

• the probabilities that the 28-day and 10-day light curves are consistent with a constant 
value. 

For the brightest GeV sources. Table [6] lists the parameters of the power-law fits to the 
low-energy (0.2 GeV-1 GeV) and high-energy (1 GeV-300 GeV) bands. For the TeV sources 
not detected by Fermi, Table [7] gives the GeV upper-limit over the Fermi energy range. 
Finally, Table [8] presents the extrapolations of the GeV spectra to the TeV domain, listing 
the differential flux extrapolated (or measured) at 100 GeV, the integral flux in the TeV 
band (0.2 TeV-10 TeV) and the photon index found by fltting the EBL-corrected spectrum 
between 100 GeV and 1 TeV with a power-law. 

4.1. TeV sources detected by the Fermi LAT 

The TeV sources detected by the Fermi LAT are discussed individually below. For 
sources with a published TeV spectrum the Fermi spectrum is compared with the TeV. 
Unless otherwise noted, when more than one TeV spectrum is available in the literature, the 
one corresponding to the lowest flux state is chosen for comparison. Since references to the 
initial detection of each TeV source and to the TeV spectra chosen are given in Tables [U and 
|2l they are not repeated in the text below. As discussed above, the TS maps, spectra and 
light curves are presented in Figures [H [21 H] and [5] and in Tables O [6] and [H and the reader 
is not directed to them individually in the discussion below. 

3C 66A/B: TeV 7-ray emission from this region was initially reported by the Crimea 
Observatory group and later, based on observations with the more sensitive VERITAS and 
MAGIC instruments, with a flux that was less than 1/100 that claimed in the original 
detection. VERITAS observations during 2007 & 2008 lead to the detection of a flare from 
3C 66A (in October 2008), an IBL, while ruling out 3C 66B, separated by 0.12° from 3C 66A, 
at a level of 4.3(T. In contrast, the MAGIC observations during 2007 are consistent with 
the emission originating from 3C 66B, a radio-galaxy at a distance oi z = 0.0211; they 
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rule out 3C 66A at a probability of 85.4%. The MAGIC observations revealed a harder 
source, with a significantly lower flux than the later VERITAS observations, which were 
taken during a flaring episode. The MAGIC observations showed no evidence for variability. 
GeV 7-ray emission from the region of 3C 66 A was discovered b y EGRET, although the 



signal was contaminated by a nearby pulsar (IKuiper et al.ll2000l ). Details of a dedicated 



multiwavelength c ampaign on 3C 66 A, involving Fermi, during the period of the VERITAS 



flare are given by iReyes et al.l (120091 ). The redshift of 3C 66A is assumed to be 2; = 0.444. 



This value, however, is based on two measurer nents of a single w eak line in the spectrum of 



the galaxy, and is considered to be uncertain (IMiller et al.lll978l ) 



The positions of 3C 66A ("x") and 3C 66B ("+") are marked in the TS map for this 
field. The emission is distributed throughout a broad region around both sources, consistent 
with the Fermi LAT PSF (which is large at low energies), however it can be seen that the 
centroid of the emission is constrained to a relatively small region containing 3C 66A (at a 
68% confidence level) and excluding 3C 66B (at >99% level). The Fermi LAT emission above 
1 GeV is well described by a power law with index of F = 1.98 ± 0.04, which is extrapolated 
to the TeV regime using the redshifts of both 3C 66A (dashed line) and 3C 66B (dot-dashed 
line). The VERITAS and MAGIC spectral measurements are also shown. It is clear that the 
spectrum extrapolated in the distant "3C 66A" scenario is in better agreement with the TeV 
measurements than the close-by "3C 66B" scenario. The latter case would require significant 
turnover in the intrinsic spectrum above 100 GeV to agree with the MAGIC measurements. 
The 28-day and 10-day lightcurves show evidence for variability, with a factor of 5-6 between 
the highest and lowest fluxes. As a result of this, and since the VERITAS measurements 
showed evidence for a flaring state, and the fact that the redshift of 3C 66A is uncertain, 
we do not claim that the extrapolated, averaged GeV spectrum is an exact match to the 
VERITAS points, only that their superficial agreement is suggestive of the dominance of a 
more distant source. Taken together the positional and spectral information indicate that 
the bulk of the GeV emission arises from 3C 66A. 

RGB J0710+591: Detected recently by the VERITAS collaboration, detailed spectral 
information at TeV energies has not yet been published for this HBL. RGB J0710+591 is 
weakly detected by the Fermi LAT, with indications of a hard spectrum. The low statistics 
at GeV energies mean that the extrapolation into the TeV regime is not constraining. 

S5 0716+714: This recently detected MAGIC source was reported with a preliminary 
flux level of 0(> 400 GeV) ~ 10~^^ cm~^s~^. Fermi-LAT observations reveal highly sig- 
nificant GeV emission, with a falling spectrum. The extrapolated GeV spectrum gives a 
flux of (l)ext{> 400 GeV) ~ 0.07 x 10~^^ cm~^s~^, indicating that the source was likely in a 
particularly bright state during the MAGIC observations. Indeed, Swift observations con- 



temporaneous with the MAGIC detect ion revealed the highest X-ray flux ever measured 
from S5 0716+714 Jciommi et aliboosh . 



lES 0806+524-' Fermi detects significant emission from this object, which is consistent 
with a flat spectrum of F = 2.04 ± 0.14. The extrapolation of this power-law to TeV 
energies agrees well with the spectrum measured by VERITAS. The VERITAS observations 
did not reveal any significant variability on timescales of months, but the flux was too low 
to probe shorter timescales. The Fermi LAT observations show only marginal evidence for 
variability on 28-day timescales, and we therefore suggest that it is reasonable to equate the 
time-averaged GeV and TeV spectra. 

lES 1011+496: The spectrum of this bright Fermi source is analyzed in two energy 
bands. The high energy band {E > 1 GeV) is consistent with a power law of index F = 
1.96 ±0.09. The TeV spectrum from MAGIC is considerably softer, with the lowest spectral 
measurement made at 150 GeV. The highest energy photon detected by Fermi has an energy 
of 168 GeV, and hence is in the range covered by the MAGIC spectrum. The flux of the two 
measured spectra are consistent in the overlapping region, and the extrapolated, absorbed 
Fermi spectrum agrees well with the measured TeV points. No evidence of variability is seen 
in the GeV or TeV observations. 

Markarian 4^1: A very bright source in the GeV regime, the spectrum of Markarian 421 
is measured with precision by the LAT. No indication of spectral curvature is found; the 
spectrum above 1 GeV is well described by a simple power law of index F = 1.78 ± 0.04. A 
similar value is found for the spectrum below 1 GeV. The highest energy photon detected 
by the LAT from Markarian 421 was reconstructed at ~ 800 GeV, and five photons were 
found with E > 150 GeV. This source has the highest degree of overlap between GeV and 
TeV spectra of any of the TeV blazars. There are a considerable number of TeV spectral 
measurements available in the literature; we adopt the MAGIC spectrum, which has the 
lowest energy threshold, was made in a relatively low flux state (O.5 0crab)- The spectrum 
during an extreme flaring state, measured by Whipple in 2001, is also shown. The differential 
flux measured by the Fermi LAT at 100 GeV is compatible with that found by MAGIC, 
nevertheless, the extrapolation of the Fermi spectrum leads to an overestimation of the 
integral flux above 200 GeV. The extrapolated photon index is 1.9, clearly harder than any 
reported in the literature. It is impossible to reconcile the GeV and TeV spectra on the 
basis of EBL absorption alone, and we conclude there is a turnover in the intrinsic spectrum 
in the neighborhood of 100 GeV, a region of falling sensitivity for both Fermi and the TeV 
instruments. 

Markarian 180: Based on detection at a level of TS = 50, the Fermi observations show 
no evidence for variability and yield a power-law spectral index of F = 1.91 ± 0.18. The 
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highest energy photon associated with Mrk 180 is 14 GeV, a decade lower than the TeV data 
points reported by MAGIC. The TeV spectrum is softer than a simple extrapolation from the 
Fermi regime, but has a larger flux at 150 GeV. MAGIC did not detect any variability from 
the object, however their observations were triggered by a particularly high optical state, 
which might indicate that the TeV spectrum is not representative of an "average" state. 

lES 1218+304: This object lies close to two others considered in this study, W Comae 
and ON 325. All three occupy a single region of interest (see section [3]) for Fermi, and must 
hence be analyzed together. In the GeV regime the brightest by far is ON 325, an LBL 
which was detected in the OFGL survey but not detected by TeV instruments. The TeV 
source lES 1218+304 lies only ~ 0.75° from ON 325, well within the PSF of the Fermi LAT, 
at least at lower energies, and W Comae, also a TeV source, lies ~ 2° away. 

Fermi detects significant emission from the region of lES 1218+304. The Fermi spec- 
trum is well described by a power law with an index of F = 1.63 ± 0.12, making it one of the 
hardest sources in the sample. In the TS map for this region, the centroid of the GeV emis- 
sion is offset by ~ 4arcmin, with the blazar located on the 95% confidence contour. Since 
the TS map shows the residual signal after the known sources have been accounted for, it 
is possible that small errors in the modeling of ON 325 (located just beyond the right edge 
of the map), possibly resulting from inaccuracies in the IRFs, introduce a systematic shift 
in the centroid of the residual emission. However, in light of the extremely hard spectrum 
measured by Fermi, we consider the GeV emission to be associated with lES 1218+304. No 
variability is detected in the Fermi light curve. During their original observations, MAGIC 
and VERITAS observed no evidence for variability at TeV energies, and the spectra they 
produced agree well. However, during observations in January and February 2009, just after 
the time period considered in this study, VER ITAS detected a fla re from lES 1218+304, dur- 



ing which its flux increased by a factor of five (jimran et al.ll2009l ). The GeV and TeV spectra 



are close to overlapping at 100 GeV, and an extrapolation from GeV energies agrees quite 
well with the TeV data points. Again, the 7-ray emission evidently peaks in the 50-150 GeV 
range. 

W Comae: The light curves on 10 and 28-day timescales show a decline of a factor of 
2-3 in flux between the start and end of the study period. A relatively bright Fermi source, 
the averaged spectrum below 1 GeV is consistent with being fiat or moderately increasing, 
while the higher energy band is softer, with F = 2.16 ± 0.10. The TeV emission has also 
proved to be highly variable, with VERITAS reporting a dramatic flare, during which their 
spectral measurements were derived. In light of the variability in the GeV and TeV regimes, 
it is not particularly surprising that the extrapolated GeV spectrum does not match well 
with the TeV data. 



-17- 



3C 279: The only FSRQ detected to date at TeV energies, 3C 279 is the strongest 
Fermi source in this study. The GeV spectrum shows clear evidence for curvature, with 
the peak of the emission lying below 200 MeV. A strong flare occurred during the period 
of this study, between MJD 54780 and 54840, during which the flux increased by a factor 
of ~ 7 and then declined. We show the spectra from the pre-flare (MJD< 54780) and 
peak-flaring (54790 <MJD< 54830) periods separately, along with the TeV measurements 
from MAGIC. The spectral indices for the E > 1 GeV components are similar in the two 
states. The predicted TeV flux from an extrapolation of the flaring GeV spectrum is 0(> 
100 GeV) = 3.5 ± 1.3 x 10^^^ cm^^s^^, an order of magnitude below the flux reported by 
MAGIC during the February 2006. The extrapolated index for both states is 3.6 ± 0.3, 
in agreement with that reported by MAGIC. In contrast, below 1 GeV the Fermi spectral 
index changed significantly between the two flux states, becoming harder (AF = 0.52 ±0.13) 
during the flaring period, suggesting that the peak of the high energy component increased 
in energy. 

PKS I424+24O: Prompted by an initial detection in the OFGL list by Fermi, this object 
was subsequently detected at TeV energies by VERITAS and confirmed by MAGIC. The 
GeV spectrum is hard (F = 1.85 ± 0.05). No redshift measurement has been made for this 
object, and, therefore, we do not extrapolate the GeV spectrum into the TeV regime. At 
the present time, no TeV spectrum has been published by VERITAS or MAGIC. 

H 1426+428: Fermi detects weak emission from the region of this HBL, with a spectral 
index of F = 1.47 ± 0.30. No evidence for variability is seen. A powerful, distant FSRQ 
(B3 1428+422, z = 4 72) lies 40 arc seconds from H 1426+42 8, too close to be resolved 



separately by Fermi (IFabian et al.lll999l : ICostamante et al.ll200ll ). However, the hard GeV 



spectral index strongly suggests that H 1426+428 is the source of the bulk of the emission 
detected by Fermi. At TeV energies, H 1426+428 was detected during active periods by 
Whipple, HEGRA, CAT and others, but no detections have been reported with the more 
sensitive third-generation lACTs. Extrapolated to TeV energies, the GeV spectrum would 
result in an integral flux of ~ 0.08</)crab, compared with 0.17(/>crab {E > 350 GeV) measured 
by Whipple and 0.08 & O.O3 0crab {E > 1 TeV) by HEGRA over two different periods. 

PG 1553+113: Detected by H.E.S.S., PC 1553+113 is one of the softest TeV sources, 
with Ftgv = 4.0 ± 0.6. In contrast, in the GeV regime, Fermi detects a bright, hard source, 
with a spectral index of F = 1.69 ± 0.04. The spectra measured by Fermi overlaps with 
those measured by H.E.S.S. and MAGIC, and are in good agreement at around 150 GeV. 
This source therefore has a strong spectral break of AF k, 2.3, which must be explained 
either through absorption with the EBL or through some mechanism intrinsic to the blazar. 



Due to the almost complete lack of measurable spectral lines, the redshift of this HBL 
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has not been established. Several indirect methods place it in the range 0.09 < z < 0.78. 
Therefore we extrapolate the flux from the GeV to TeV energy ranges by assuming both the 
lower and upper limit on redshift. In the case oi z = 0.78, the EBL absorption is sufficient 
that the extrapolated GeV spectrum is in good agreement with the TeV measurements, 
leading credence to the hypothesis that the redshift of PC 1553+113 is significantly larger 
than the lower bound of the allowable redshift range. Indeed, PG 1553+113 might be the 
most distant TeV object detected to date. This source is the subject of an independent 



Fermz-LAT paper flAbdo et al.ll2009gf ). 



Markarian 501: Many episodes of fiaring have been detected from Mrk 501 with previous 
generation TeV instruments. To date, however, the spectrum for a low-level fiux state has 
not been published by third-generation instruments. In the GeV regime, EGRET detected 
emission from Mrk 501 after the initial discovery at TeV energies, and significant emission 
is observed by the LAT. The GeV spectrum of Mrk 501 is well fitted with a simple power 
law with F = 1.79 ± 0.06. In contrast to the historical TeV emission, the GeV fiux shows 
no evidence of variability. With Fext = 1-86 ± 0.05, the extrapolated GeV spectrum is 
harder than all spectral measurements made with TeV instruments, indicating the presence 
of curvature in the intrinsic spectrum of the source. TeV instruments have long detected 
evidence for curvature in the spectrum of this object at 1 TeV. The GeV spectrum is shown 
with TeV measurements made during a moderate fiare (O.40crab, E > 150 GeV) and during 
a very high state (1.8 0crab, E > 500 GeV). 

lES 1959+650: Fermi detects emission from this object with a fiat spectrum, finding no 
evidence of variability over the period of this study. At TeV energies, lES 1959+650 had long 
been detected only while fiaring. However, during a dedicated multiwavelength campaign 
in 2006, MAGIC measured its spectrum in a low fiux state. An extrapolation of the GeV 
butterfiy overlaps these MAGIC measurements indicating that the underlying spectrum is 
largely compatible with a single power law over the full 7-ray regime. The difference between 
the measured and extrapolated spectral indices, however, is Fxev — Fext = 0.41 ± 0.19, 
indicating that there is some evidence for curvature between the two bands (at the 2a level). 

PKS 2005-489: A southern hemisphere HBL detected in the TeV domain by H.E.S.S., 
this source is one of the softer TeV blazars with Fxev = 4.0 ± 0.4. Fermi detects significant, 
hard GeV emission with an index of F = 1.91 ± 0.09. No evidence of variability is seen 
by Fermi on timescales of months while H.E.S.S. observes variability only on timescales 
longer than a year. The difference between the indices of the extrapolated GeV and H.E.S.S. 
spectra is AF = 1.8 ± 0.4, indicating a clear break at a few hundred GeV. 

PKS 2155-304: The results of a dedicated multiwavelength campaign on this object, 
including GeV and TeV observations with Fermi and H.E.S.S. and simple SSC modeling. 
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are reported by lAharonian et al.l (l2009d ). Since the publication of these results, an improved 
set of IRFs have become available, which correct for an overestimate of the effective area at 
lower energies that results from on-orbit "pile-up" of events in the detector. A reanalysis 
of the Fermi data with these IRFs results in an increase in flux of ~ 15% with almost no 
change in the spectral index. During the 5.5 month period of this study, the GeV flux of 
PKS 2155-304 varied by a factor of ~ 3, with a maximum of 145 x 10~^cm~^s~^, ~ 1.5 
times higher than the average flux. The change in spectral index between the GeV and TeV 
measurements can only partly be explained by EBL absorption, the remainder presumably 
resulting from some process intrinsic to the source. 

BL Lacertae: The spectral index measured by Fermi from this LBL is relatively soft, 
with r = 2.43 ±0.1. No evidence for variability is seen over the period of the study. The 
extrapolated flux is approximately one third that measured by MAGIC during a flaring 
episode. Since this flux represents an estimate of the TeV flux in the optimistic case that 
there is no intrinsic curvature, we conclude that the low flux state from this object will likely 
be difficult to measure at TeV energies without a significant investment of observing time. 

lES 2344+514: One of the fainter sources in this study, lES 2344+514 is detected 
with a TS of only 37 and the highest energy photon collected has E = 53 GeV. The Fermi 
spectrum is not consistent with the MAGIC spectrum obtained in a low state, but the poor 
statistics are insufficient to make a reliable prediction of the TeV flux. 

M 87 & Centaurus A: These nearby sources are two of the three radio galaxies thus 
far detected at TeV energies. Both are classifled as FRl, thought to be the parent class 
of BL Lacs. While Cen A has not shown evidence of variability on time scales of days or 
months, M 87 has undergone several flaring episodes on time scales as short as a day. In the 
GeV domain, Fermi detects both of these sources. The flux of M 87 between 0.2 GeV and 
300 GeV, is ten times lower than the flux of Cen A and is too faint to make strong predictions 
for the TeV emission. In contrast, the detection of Cen A with TS = 308 yields a strongly 
constraining extrapolation to the TeV domain. The Fermi butterfly underestimates the TeV 
measurements from Cen A by a factor 10. M 87 and Cen A are the subject of dedicated 



Fermi papers flAbdo et al.ll2009(] j3l 



4.2. TeV sources not detected by the Fermi LAT 

Of the 28 TeV-detected AGN, only the six listed in Table [7] (and 3C 66B, for which 
we cannot produce upper limits due to contamination by 3C 66A) have not been detected 
at GeV energies in 5.5 months of data taking with Fermi. The differential fluxes of these 
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objects at 200 GeV, calculated from measured TeV spectra (column F200 in Table [2]), are all 
below the median flux, med(F2oo) = 0.171, from the full sample. In fact, they are six of the 
eight TeV objects with the smallest F200 fluxes. 

Upper limits are calculated from the GeV observations assuming two scenarios for the 
spectral index, Fhard = 1-5 and Fsoft = 2.0, and are given in Table [71 These hmits are 
extrapolated into the TeV regime in the usual way, and are presented in Figure [31 In 
general, these two sets of extrapolated limits bracket the TeV measurements. 



4.3. GeV sources with upper limits in the TeV regime 

In addition to TeV-detected sources, Fermi detects emission from 17 AGN for which 
only upper limits exist at TeV energies. These GeV detections can be extrapolated to predict 
the flux that might be observable by TeV observatories, with the caveat that assuming the 
intrinsic spectrum is described by a single power-law up to the TeV regime can lead to overly 
optimistic predictions for detection. For most of the TeV sources, the spectra in the TeV 
regime are well described by a power law with an index no harder than F = 2, and it is 
reasonable to think that this will be the case for future detections. In this context, sources 
with a predicted photon index harder than F = 2 should exhibit an intrinsic spectral break. 

lES 0033+595: Observed by Whipple for 12 hours without detection, this HBL, visible 
only from the northern hem isphere, has an uncertain re dshift. Perlman tentatively measured 



0.086 (see comment in 



Falomo fc KotilainenI Eoool ) . while others have argued that it is 



more distant (e.g. ISbarufatti et al.l [20051 . who claim z > 0.24). We adopt the tentative, 
direct measurement. The Fermi spectrum is consistent with a flat power law with F = 
2.00 ± 0.13) and photons detected up to 150 GeV, resulting in an extrapolated index of 
Fext = 2.35 ± 0.08. The Whipple upper limit, 0TeV <O.166 0crab, is in agreement with 
the prediction of 0.053 (Acrah obta ined here. To reach this flux level with the sensitivity of 



VERITAS (see [Holder et al.[[2008[ ) would require an observation of approximately 2.5 hours. 



PKS B0521-365, PKS 0829+046 & 3C 273: The objects have an extrapolated flux less 
than O.O10crab- In the absence of signiflcant flaring, it will likely be difficult to detect TeV 
emission from them unless a large amount of time is dedicated to their observation. 

lES O647+25O: The Fermi LAT detects the emission of this object with a very hard 
spectral index of F = 1.66 ± 0.15 and very high maximum photon energy of 257GeV. The 
source is visible in the northern hemisphere and was observed by HEGRA for 4 hours. The 
predicted flux from the Fermi data is 0.0543 ^crab, consistent with the HEGRA upper limit 
and corresponding to 2.5 hours of observations with VERITAS. 
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lES 1028+511, lES III8+424 & IZw 187: These AGN lie outside the 95% probabihty 
contour for the origin of the emission on the TS maps. Given that there are 17 AGN detected 
in this category, the chance probabihty of discovering three outside the 95% contours is 
-PBm(> 3, 17, 1 — 0.95) = 0.050, equivalent to a 2. Oct Gaussian event. The possibility that the 
emission detected from one or more of these three regions is not associated with the particular 
AGN under study cannot be discounted. For completeness, the chance probability based on 
the 38 AGN detected in this study is PBin(> 3, 38, 1 — 0.95) = 0.30, not unreasonable, but 
begging the question as to why all such objects belong to this category (TeV non-detected 
sources in Table [5]). 

lES 1440+122: This Fermi-LAT source is formally the hardest detected in this study. 
However, since the LAT detected only 10 photons with energy E > 1 GeV, more observations 
are required before firm statements can be made about the spectral index. With a redshift of 
z = 0.162, the extrapolated index between 200 GeV and 1 TeV is 1.68, which is harder than 
any TeV source yet detected. The predicted flux of 0.231 0crab is likely overestimated by a 
large factor and a turnover, not accounted for by EBL absorption only, is required to explain 
the discrepancy between the predicted flux and the H.E.S.S. upper limits of O.O3 0crab- 



5. Evolution of detected GeV— TeV spectra with redshift 

In the LBAS study, the dependence of GeV spe ctral index of the FSRQ and BL Lac 



populations on redshift was presented (Figure 11 of lAbdo et al.l l2009al ). Although it was 
found that the two populations have different spectral properties, no significant relation 
between the gamma-ray photon index and redshift was found within each source class. The 
GeV-TeV sources provide a population in which the effects of spectral evolution with redshift 
can be studied across a much wider energy range than LBAS (although, admittedly in a 
much smaller redshift range). The presence of a redshift-dependent spectral break in these 
sources could be indicative of the effects of absorption on the EBL, and provide experimental 
evidence for this absorption in a manner independent of any specific EBL-density model. 
Appendix |A] discusses the relationship between EBL absorption and the redshift dependent 
change in spectral index between the GeV and TeV ranges, 

dT{E,z 



AF = FxeV - FceV > FtgY - Fint = S{z, E* 



(2) 

E=E* 



dlogE 

i.e. in the presence of EBL absorption, and assuming that the intrinsic spectra do not get 
harder with energy, measured spectral breaks must lie in a region of the (AF, z) plane defined 
by AF > 6{z, E*), with 6{z, E*) defined by the EBL optical depth. 

Figure [6] depicts the difference in the measured TeV and GeV spectral indices, AF, as 
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a function of the redshift, z, for 15 of the 21 extragalactic GeV-TeV source^. It is evident 
that the difference between the GeV and TeV spectral indices increases with redshift. At 
low redshifts the radio galaxies M 87 and Cen A have AF ^ 0, as do the near-by BL Lacs. 
At redshifts greater than 0.1, all of the BL Lacs are consistent with AF > 1.5. 

Pearson's correlation factor, r, is widely used to quantify the correlation between two 
variables. However, this correlation factor does not take measurement errors into account. 
Thus, to evaluate the significance of the correlation, a series of Gaussian random variables 
centered on the values of the measured spectral changes, AFj, and with width equal to their 
measurement errors are used to generate multiple simulated data sets. The full width at 
half maximum of the distribution obtained for r gives an estimate of the erroi0. The value 
obtained is r = 0.76 ± 0.14, which indicates a clear correlation. To check the robustness of 
this result, the Kendall rank, defined by 

2S 

tk 



N -{N -1)'' 



where N is the size of the data set (A^ = 15), is been calculated. For each pair of points 
from the dataset (AF,2;), those in the same order are assigned a value of +1, and the 
others assigned —1. T he sum, S, of these A^ • (A^ — l)/2 combinations is then constructed 



( iGleissner et al.ll2004l . and references therein), giving a value of tk, ranging from — 1 < tj^ < 
+ 1, indicating the degree of correlation. The error on tk is calculated in the same way as 
for r. The same conclusion is established with this test, tk = 0.68 ± 0.15. 

The effects of systematics present in Figure [6] must be considered when determining the 
validity of concluding that it reflects the effects of EBL absorption. Such systematics are 
difficult to evaluate in a quantitative manner, appendix [B] discusses possible contributions 
in more detail. In light of these difficulties it cannot be claimed with 100% certainty that 
the observed deficit of sources at AF ^ for large redshift is a real effect. We anticipate 
that the TeV AGN not detected by the LAT in this study will ultimately be detected and 
that this figure will then be limited only by the selection bias at TeV energies. Finally, 
EBL absorption is not the only effect that could cause a redshift dependent spectral break, 
evolution of the intrinsic spectra with redshift also cannot be excluded, although this may be 
difficult to reconcile with the LBAS study, which did not detect any evolution in the larger 
LAT BL Lac sample. 



^We exclude 3C 66A and PG 1553+113 as z is not known, RGB J0710+59, S5 0716+714 and 
PKS 1424+240 as FxeV has not been published and 3C 279 as the GeV and TeV states are badly mis- 
matched. 

^The analytic expression for the error gives r — 0.76 ± 0.16, consistent with the Monte Carlo method. 
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6. Conclusions 

In 5.5 months of observation the Fermi LAT has detected GeV emission from 21 TeV- 
detected AGN, and from 17 AGN previously observed by TeV groups and for which upper 
hmits have been pubhshed at TeV energies. Whereas EGRET detected only a small number 
of such TeV sources, and detected those only with integration times of months and years, 
Fermi has detected the majority of the TeV blazars in its first few months of operation. Fermi 
observations will help TeV observatories optimize the limited observation time available each 
year, and will be uscftil in evaluating possible targets for observation. Fermi has sufficient 
sensitivity to participate meaningfully in simultaneous multi-wavelength campaigns and to 
measure spectra and light curves from the brighter blazars with a resolution of days to 
months. In future campaigns, it can be expected that the high-energy emission will be as 
well covered by Fermi and the TeV observatories as the low-energy peak is by instruments 
in the radio to X-ray bands. 

Many of the TeV sources exhibit an increasing spectrum (F < 2) in the GcV range 
confirming the presence of a high energy peak in z/Fj, representation. This is the first large- 
scale characterization of the full 7-ray emission component for this class of energetic AGN. 
More detailed modeling of these 7-ray sources, which is beyond the scope of this paper, will 
become possible as more data are acquired by Fermi and the flight-cahbrated IRFs become 
available (extending the effective energy below 100 MeV). The MAGIC-II and H.E.S.S. II 
instruments will increase the range over which the sensitives of TeV instruments overlap 
with Fermi at lower energies, producing better coverage at energies between 50 and 200 GeV 
where a number of TeV sources seem to have turnovers in their measured spectra. 

The intrinsic spectrum for some of the TeV sources can be well described by a single 
power-law across the energy range spanned by the Fermi LAT and the TeV observatories, 
with any breaks in the measured 7-ray spectra between the two regimes being consistent 
with the effects of absorption with a model of minimal EBL density. For other objects, 
however, a softening of the intrinsic spectrum is required to match the TeV measurements. 
This could be due to softening intrinsic to the IC component - itself reflecting curvature in 
the relativistic electron or seed photon distributions. 

Based on an extrapolation of the GeV spectra measured by Fermi to TeV energies, a 
number of previously observed TeV candidates are good candidates for re-observation with 
TeV instruments: lES 0033+595 and lES 0647+250 are two such targets. 

Redshift-dependent evolution is detected in the spectra of objects detected at GeV and 
TeV energies. The most reasonable explanation for this is absorption on the EBL, and as 
such, it would represent the flrst model-independent evidence for absorption of 7 rays on 
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the EBL. Future observations with Fermi and TeV instruments have the potential to probe 
t{E,z) in a more quantitative manner. 
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A. Relationship between GeV— TeV break index and EBL 



Writing the measured TeV spectrum as FTev(-E'), the unabsorbed (intrinsic) spectrum 
as Fint(-E) and the redshift-dependent optical depth due to EBL absorption as t{E, z), the 
effects of the EBL on the measured spectrum is given by. 



(Al) 



If both the measured and intrinsic spectra can be approximated as power laws over the 
range of the TeV observations, FTev{E) ^ CTev(^/^o)"^^"^ and Fi^t{E) ^ Cint(^/£^o 
then the EBL effects must also be given by a power law, say e~'^^^'^^ ~ Cr{E /Eo)~^^^'^*\ 
The power-law index of the EBL absorption, 6{z,E*), is a function of the redshift and the 
energy range over which the TeV observations are made (denoted for convenience as a de- 
pendence on some energy, E*, at which the measured TeV data most constrains the fitted 
spectrum). Equating the two express ions for the absorption to first order in x = log{E/E*) 



gives 6{z,E*) ^ dr^x, z)/dx\x=o (see IVassilievI I2OOOI : IStecker fc ScuUyl l2006l . for further dis- 
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cussion of linear and polynomial expansions of the EBL optical depth). Equation lAll relates 
6{z, E*) to the intrinsic and measured spectral indices, 

CTME/Eor^-'^- = CrCUE/E,r'^^'''*^-^^-\ or 

5{z, E*) = TTev - Tint ~ Ttcv - TceV, (A2) 

where it is further assumed that the measured GeV spectral index can be used as a proxy 
for the intrinsic index in the TeV regime. This is equivalent to assuming that there is no 
curvature in the intrinsic spectrum between the GeV and TeV energy regimes and that 
absorption on the EBL does not affect the spectrum in the LAT energy range, which is true 
for all reasonable EBL models for sources with z < 0.5. For GeV-TeV detected sources, 
this suggests that the variable AF = Ftev — Fqcv probes the effects of EBL absorption in a 
manner independent of any specific model of EBL density. 

For real sources there is curvature in the intrinsic spectrum between the GeV and TeV 
regimes, and the measured Fqev is not a perfect estimator for Fint in the TeV regime. In 
general it has been found that the curvature in the differential spectra is concave, and the 
intrinsic spectrum at TeV energies is expected to be softer than the measured GeV spectrum. 
Therefore it is expected that for real sources: 

AF = FTeV - FceV > FxeV " F^^ = 5{z, E*), (A3) 

so that sources would occupy a region of the space of (AF, z) above some curve defined by 
the EBL. 



B. Discussion of systematics in evolution of spectra with redshift 

A number of components contribute to the systematics in Figure E] and are addressed 
below: (i) systematic errors on the points themselves and (ii) the effects of the criteria used 
to select targets for the study. 

The systematic errors on the measurements of Fcev for the individual sources are pre- 
sented in Table El and are generally smaller than statistical errors, with the largest being 
~ 0.25. Similarly, the TeV groups estimate an d report systema t ic erro rs on Fxev for each 



of the detected TeV sources. See, for example, lAharonian et al.l (j2006d ) for a discussion of 
systematic error estimation with H.E.S.S. In general, the GeV and TeV systematic errors 
are too small to explain the trend in Figure El in which AF changes by > 2.0 over the range 
of redshift in question. 
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The sources in this study are subject to a two-stage selection process, which may lead to 
regions of the phase-space of (AF, z) being inaccessible due to limitations imposed by of the 
sensitivities of Fermi and the TeV instruments. This could, in turn, lead to a false correlation 



being evident in the data - for example, see Figure 7 from the LBAS study ( lAbdo et al. 



2009al ) for a correlation which might incorrectly be claimed based on a sensitivity-limited 



sample. Targets were originally selected by TeV astronomers for observation with TeV 
instruments. Given the sensitivities of those instruments and the amount of time dedicated to 
each target, some fraction were detected, leading to 28 sources used in this study. The Fermi 
sensitivity further restricts the sample to the 15 GeV-TeV sources displayed in Figure [61 
Given the complexity (and randomness) of the selection process, it is almost impossible to 
quantify where its "sensitivity" limits are in the space of (AF, 2;). In the null hypothesis, 
that there is no EBL effect present in the results of Figure El we attempt to evaluate in a 
qualitative manner whether the lack of sources at AF ^ for larger redshifts could arise from 
a selection effect. The primary selection of targets is based on detection at TeV energies. 
Sources in this region of the plot would have harder TeV spectra than those actually detected 
at the larger redshifts. Since, it is unlikely that TeV astronomers are deliberately biasing 
the sample toward softer distant AGN, the major effect producing this bias would have to 
result from the sensitivity limit of the instrument. The instrumental sensitivity directly 
limits the space of (0, F), requiring > 4>um(r) for detection (as in the LBAS figure cited 
above), and these limits transfer to the space of (AF, z) only through convolution with the 
source function /(0,F;z) = d^P/d0dF. Therefore, a sharp cut-off in this space should 
not be expected, rather a slower decrease in source counts into the "forbidden" region. In 
the no-EBL hypothesis, and further assuming that evolution in the source function, the 
primary consideration is whether the decrease in measured flux with distance coupled with 
the sensitivity limit would lead to an evolution in the population of detectable sources with 
redshift. There is very little published material addressing the flux sensitivity of current TeV 
instruments as a function of spectral index, however it is reasonable to presume that TeV 
instruments are more sensitive to sources with harder spectra, than to those with softer: 
they have better background rejection and an improved PSF at higher energies. Indeed, 
TeV instruments often use "hard cuts" to improve the sensitivity for hard sources. In this 
case, it would be expected that TeV instruments would preferentially detect harder sources 
(with lower fluxes) at larger redshifts, whereas this is exactly the opposite of what is actually 
observed, i.e. that the majority of the distant AGN are softer - no hard, weak TeV AGN 
has been detected to date. 
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Table 1. AGN detected at TeV energies. 



Name 


aj2000 


<5.J2000 


lype 




z 


rCCr 


Blazars: 






















RGB J0152+017 


01'' 


52'" 


39.6= 


+01° 


47' 


17" 


HBL 





080 


1 


3C 66A 


02'' 


22"" 


39.6= 


+43° 


02' 


08" 


IBL 





444b 


2,3"= 


lES 0229+200 


02'' 


32"" 


48.6= 


+20° 


17' 


17" 


HBL 





140 


4 


lES 0347-121 


03'' 


49"' 


23.2= 


-11° 


59' 


27" 


HBL 





188 


5 


PKS 0548-322 


05'' 


50"' 


40.8= 


-32° 


16' 


18" 


HBL 





069 


6 


RGB J0710+591 


07'' 


10™ 


30.1= 


+59° 


08' 


20" 


HBL 





125 


7 


S5 0716-1-714 


07'' 


21™ 


53.4= 


+71° 


20' 


36" 


LBL 





300 


8 


lES 0806-1-524 


08'' 


09™ 


49.2= 


+52° 


18' 


58" 


HBL 





138 


9 


lES 1011-1-496 


10'' 


15™ 


04.1= 


+49° 


26' 


01" 


HBL 





212 


10 


lES 1101-232 


11'' 


03'" 


37.6= 


-23° 


29' 


30" 


HBL 





186 


11 


Markarian 421 


11'' 


04™ 


27.3= 


+38° 


12' 


32" 


HBL 





031 


12 


Markarian 180 


11'' 


36™ 


26.4= 


+70° 


09' 


27" 


HBL 





046 


13 


lES 1218-1-304 


12'' 


21™ 


21.9= 


+30° 


10' 


37" 


HBL 





182 


14 


W Comae 


12'' 


21™ 


31.7= 


+28° 


13' 


59" 


IBL 





102 


15 


3C 279 


12'' 


56™ 


11.2= 


-05° 


47' 


22" 


FSRQ 





536 


16 


PKS 1424+240 


14h 


27™ 


00.4= 


+23° 


48' 


00" 


IBL 






17 


H 1426+428 


14'' 


28™ 


32.7= 


+42° 


40' 


21" 


HBL 





129 


18 


PG 1553+113 


15'' 


55™ 


43.0= 


+11° 


11' 


24" 


HBL 





09 - 0.78 


19 


Markarian 501 


16'' 


53™ 


52.2= 


+39° 


45' 


37" 


HBL 





034 


20 


lES 1959+650 


19'' 


59™ 


59.9= 


+65° 


08' 


55" 


HBL 





048 


21 


PKS 2005-489 


20'' 


09™ 


25.4= 


-48° 


49' 


54" 


HBL 





071 


22 


PKS 2155-304 


21'' 


58™ 


52.1= 


-30° 


13' 


32" 


HBL 





117 


23 


BL Lacertae 


22'' 


02™ 


43.3= 


+42° 


16' 


40" 


LBL 





069 


24,25'= 


lES 2344+514 


23'' 


47™ 


04.8= 


+51° 


42' 


18" 


HBL 





044 


26 


H 2356-309 


23'' 


59™ 


07.9= 


-30° 


37' 


41" 


HBL 





167 


27 


Others 






















3C 66B 


02'' 


23™ 


11.4= 


+42° 


59' 


31" 


FRl 





02106 


28 


M 87 


12'' 


30™ 


49.4= 


+12° 


23' 


28" 


FRl 





004233 


29 


Centaurus A 


13'' 


25™ 


27.6= 


-43° 


01' 


09" 


FRl 





00183 


30 



^See notes for Table |3] for explanation of object types. 
''The redshift of 3C 66A is considered to be uncertain. 

'^Detection of E>lTeV em ission from 30 66A and BL Lacertae was first claimed 
bv lNeshpor et al ] ||l998l . l200ih . The measured fluxes are not consistent with the later 
measurements made with more sensitive instruments. 
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Table 2. Flux (0), photon index (F) from measurements of AGN with TeV instruments, 
along with threshold energy (-Ethres) of the observation. The differential flux at 200 GeV 
(-P200) is also calculated from the TeV spectrum for comparison between objects. See text 

for further details. 



Name 


-^thrcs 


0{> -Ethrcs) 


r 


-F200 


Note 


Ref 




[GcV] 


[10""cm-2s-i] 


[1] 


[lO-'-'cm-^s-iTeV-i] 






Blazars: 














RGB J0152+017 


300 


0.27 ± 0.05 


2.95 + 0.36 


0.058 


No variability 


1 


3C 66A 


200 


1.3 ± 0.1 


4.1 + 0.4 


0.201 


Flaring state 


2 


lES 0229+200 


580 


0.094 ± 0.015 


2.50 + 0.19 


0.003 


No variability 


3 


lES 0347-121 


250 


0.39 ± 0.01 


3.10 + 0.23 


0.065 


No variability 


4 


PKS 0548-322 


200 


0.33 ± 0.07 


2.8 + 0.3 


0.030 




5 


RGB J0710+591 


300 


RJ 0.016 0Crab 






ATEL 1941 


6 


S5 0716-1-714 


400 


~ 1 <ACrab 






ATEL 1502 


7 


lES 0806-1-524 


300 


0.22 ± 0.05 


3.6 + 1.0 


2.231 


Low flux state? 


8 


lES 1011-1-496 


200 


1.58 ± 0.32 


4.0 + 0.5 


0.237 


Flaring state 


9 


lES 1101-232 


225 


0.52 ± 0.14 


2.94 + 0.20 


0.063 


No variability 


10 


Markarian 421 


383 


573.3 ± 57.9 


2.31 + 0.04 


87.96 


Flaring state 


11 




200 


26.2 ± 2.1 


2.20 + 0.08 


1.572 


Lowest flux state 


12 


Markarian 180 


200 


2.3 ± 0.7 


3.3 + 0.7 


0.264 


Flaring state 


13 


lES 1218-1-304 


200 


1.22 ± 0.26 


3.08 + 0.34 


0.125 




14 


W Comae 


200 


1.99 ± 0.07 


3.81 + 0.35 


0.280 


Flaring state 


15 


30 279 


100 


51.5 ± 8.2 


4.11 + 0.68 


0.931 


Flaring state 


16 


PKS 1424+240 


200 


Si 0.02 0Crab 






ATEL 2084 


17 


H 1426+428 


280 


2.4 + 4.1 


3.50 + 0.35 


0.696 




18,19^ 


PG 1553+113 


200 


4.8 + 1.3 


4.0 + 0.6 


0.720 


Low flux state? 


20 


Markarian 501 


300 


30.3 ± 1.9 


2.22 + 0.04 


3.03 


Flaring state 


21 




150 


12.4 + 0.8 


2.45 + 0.07 


0.592 


Lowest flux state 


22 


lES 1959+650 


1300 


2.50 + 0.46 


2.83 + 0.14 


7.030 


Orphan flare 


23 




150 


3.42 + 0.92 


2.58 + 0.18 


0.171 


Low flux state 


24 


PKS 2005-489 


200 


0.62 + 0.1 


4.0 + 0.4 


0.093 




25 


PKS 2155-304 


200 


172 + 0.05 


3.19 + 0.02'= 


20.6 


Flaring state 


26 




300 


4.2 + 0.75 


3.32 + 0.06 


1.248 


Low flux state 


27 


BL Lacertae 


200 


0.6 + 0.2 


3.6 + 0.5 


0.078 




28 


lES 2344+514 


200 


2.39 + 0.3 


2.95 + 0.12 


0.233 


Low flux state 


29 


H 2356-309 


200 


0.41 + 0.05 


3.09 + 0.24 


0.043 


No variability 


30 


Others 














3C 66B 


150 


0.73 + 0.15 


3.1 + 0.31 


0.042 


Flaring state 


31 


M 87 


730 


0.025 + 0.03 


2.62 + 0.35 


0.016 




32 


Centaurus A 


250 


0.156 + 0.067 


2.73 + 0.45 


0.019 


Low flux state? 


33 



''The flux reported by the HEGRA collaboration is 0.08 0Crab above 1 TeV. 

''A broken power law with Pi = 2.71 + 0.06, r2 = 3.53 + 0.05 and Bbroak = 430 + 22 GeV was preferred. 
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Table 3. AGN with published upper hmits at TeV energies. 













TeV limit 




Name 


"J2000 


'5j2000 


Type^ 


z 


Flux 

1 


Energy 
[GcV J 


Ref 


III Zw 2 


OO'' 10™ 31.0^ 


+10° 58' 30" 


Syl 


0.089 


<0.027 


>430 


1 


lES 0033+595 


00'' 35™ 52.6^ 


+59° 50' 05" 


HBL 


0.086'' 


<0.11 


>390 


2 


NGC 315 


OQh 57m 48 gs 


+30° 21' 09" 


FRl 


0.016 


<0.05 


>860 


3 


4C+31.04 


01'' 19™ 35.0= 


+32° 10' 50" 


RG 


0.060 


<0.14 


>760 


3 


lES 0120+340 


01'' 23™ 08.6= 


+34° 20' 49" 


HBL 


0.272 


<0.032 


>190 


4 


lES 0145+138 


01'' 48™ 29.8= 


+14° 02' 19" 


HBL 


0.125 


<0.015 


>310 


5 


UGC 01651 


Q2h ggm 33 Qs 


+35° 47' 50" 


RG 


0.038 


<0.07 


>790 


3 


BWE 0210+1159 


02'' 13™ 05.2= 


+12° 13' 11" 


LBL 


0.250 


<0.012 


>530 


1 


RGB J0214+517 


02'' 14™ 17.9= 


+51° 44' 52" 


HBL 


0.049 


<0.17 


>430 


6 


PKS 0219-164 


02'' 22™ 01.0= 


-16° 15' 17" 


LBL 


0.698 


<0.27 


>1780 


3 


NGC 1054 


Q2h 42m J5 ys 


+18° 13' 02" 


Sy 


0.033 


<0.02 


>860 


3 


NGC 1068 


Q2h 42m 40.7s 


-00° 00' 48" 


Sy2 


0.004 


<0.013 


>210 


5 


V Zw 331 


03'' 13™ 57.6= 


+41° 15' 24" 


LBL 


0.029 


<0.09 


>870 


3 


NGC 1275 


03'' 19™ 48.2= 


+41° 30' 42" 


FRl 


0.018 


<0.03 


>850 


3 


RX J0319.8+1845 




+18° 45' 34" 


HBL 


0.190 


<0.033 


>190 


4 


B2 0321+33B 


03h 24m 2= 


+34° 10' 46" 


NLSyl 


0.063 


<0.10 


>400 


7 


lES 0323+022 


03'' 26™ 13.9= 


+02° 25' 15" 


HBL 


0.147 


<0.015 


>210 


5 


4C +37.11 


04h 05m 4g.3s 


+38° 03' 32" 


RG 


0.055 


<0.05 


>800 


3 


lES 0414+009 


04h iQui 52.3s 


+01° 05' 54" 


HBL 


0.287 


<0.057 


>230 


4 


3C 120 


04'' 33"' 11.1= 


+05° 21' 16" 


RG 


0.033 


<0.004 


>230 


5 


MG J0509+0541 


05'' 09™ 26.0= 


+05° 41' 35" 


IBL 




<0.11 


>960 


3 


4C+01.13 


05'' 13™ 52.5= 


+01° 57' 10" 


BL Lac 


0.084 


<0.10 


>1010 


3 


Pictor A 


05'' 19™ 49.7= 


-45° 46' 45" 


FR2 


0.034 


<0.014 


>220 


5 


PKS B0521-365 


05'' 22™ 58.0= 


-36° 27' 31" 


LBL 


0.055 


<0.042 


>310 


1 


EXO 0556.4-3838 


05'' 58™ 06.2= 


-38° 38' 27" 


HBL 


0.302 


<0.051 


>220 


5 


PKS 0558-504 


05'' 59™ 47.4= 


-50° 26' 52" 


NLSyl 


0.137 


<0.018 


>310 


1 


lES 0647+250 


06'' 50™ 46.6= 


+25° 03' 00" 


HBL 


0.203 


<0.13 


>780 


3 


UGC 03927 


07'' 37™ 30.1= 


+59° 41' 03" 


UnC 


0.041 


<0.09 


>1090 


3 


3C 192.0 


08'' 05™ 35.0= 


+24° 09' 50" 


RG 


0.060 


<0.20 


>930 


3 


RGB J0812+026 


08'' 12™ 01.9= 


+02° 37' 33" 


BL hax: 




<0.031 


>220 


5 


3C 197.1 


08'' 21™ 33.6= 


+47° 02' 37" 


FR2 


0.130 


<0.05 


>960 


3 


PKS 0829+046 


08'' 31™ 48.9= 


+04° 29' 39" 


LBL 


0.174 


<0.06 


>1000 


3 


NGC 2622 


08'' 38™ 11.0= 


+24° 53' 43" 


Syl 


0.028 


<0.05 


>400 


7 


lES 0927+500 


09'' 30™ 37.6= 


+49° 50' 26" 


HBL 


0.188 


<0.052 


>230 


4 


84 0954+65 


09'' 58™ 47.2= 


+65° 33' 55" 


LBL 


0.368 


<0.096 


>300 


6 


MS1019.0+5139 


10'' 22™ 12.6= 


+51° 24' 00" 


BL hax: 


0.141 


<0.07 


>920 


3 


lES 1028+511 


10'' 31™ 18.5= 


+50° 53' 36" 


HBL 


0.361 


<0.29 


>400 


6 


RGB J1117+202 


11'' 17™ 06.3= 


+20° 14' 07" 


HBL 


0.139 


<0.030 


>610 


5 


lES 1118+424 


11'' 20™ 48.0= 


+42° 12' 12" 


HBL 


0.124 


<0.12 


> 430-500 


6 


Markarian 40 


11'' 25™ 36.2= 


+54° 22' 57" 


Syl 


0.021 


<0.21 


>430 


6 


NGC 3783 


11'' 39™ 01.7= 


-37° 44' 19" 


Syl 


0.010 


<0.025 


>220 


5 


NGC 4151 


12'' 10™ 32.6= 


+39° 24' 21" 


Syl.5 


0.003 


<0.07 


>790 


3 


lES 1212+078 


12'' 15™ 11.2= 


+07° 32' 05" 


HBL 


0.136 


<0.17 


>920 


3 


ON 325 


12'' 17™ 52.1= 


+30° 07' 01" 


LBL 


0.130 


<0.22 


>400-430 


6 
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Table 3 — Continued 



TeV limit 



Name 


JzUUU 




TvDe^ 


z 


Flux 


\\J^ V J 


Ref 


3C 273 


121^ 29™ 06.7= 


+02° 03' 09" 


FSRQ 


0.158 


<0.014 


>300 


1 


MS 1229.2+6430 


121^ 31™ 31.4= 


+64° 14' 18" 


HBL 


0.164 


<0.17 


> 300-430 


6 


lES 1239+069 


121^ 41™ 48.3= 


+06° 36' 01" 


HBL 


0.150 


<0.20 


> 400-430 


6 


lES 1255+244 


12^57^31.9= 


+24° 12' 40" 


HBL 


0.141 


<0.11 


>350-500 


6 


RGB J1413+436 


14^ 13™ 43.7= 


+43° 39' 45" 


RG 


0.089 


<0.06 


>400 


7 


RX J1417.9+2543 


141^17™ 56.7= 


+25° 43' 26" 


HBL 


0.237 


<0.023 


>190 


4 


OQ530 


14h 19™ 46.6= 


+54° 23' 15" 


LBL 


0.151 


<0.058 


>300 


6 


lES 1440+122 


14^ 42™ 48.3= 


+12° 00' 40" 


HBL 


0.162 


<0.033 


>290 


5 


RGB J1629+401 


29™ 01.3= 


+40° 08' 00" 


NLSyl 


0.271 


<0.09 


>400 


7 


RX J1725.0+1152 


17!^ 25™ 04.4= 


+11° 52' 15" 


HBL 


0.018 


<0.046 


>190 


4 


I Zw 187 


17*1 28™ 18.6= 


+50° 13' 10" 


HBL 


0.055 


<0.086 


> 300-350 


6 


lES 1741+196 


17!^ 43™ 57.8= 


+19° 35' 09" 


HBL 


0.083 


<0.053 


>350-500 


6 


3C 371 


181^ 06™ 50.7= 


+69° 49' 28" 


LBL 


0.051 


<0.19 


>300 


6 


Cyg A 


19!^ 59™ 28.4= 


+40° 44' 02" 


FR2 


0.056 


<0.03 


>910 


3 


PKS 2201+04 


04™ 17.7= 


+04° 40' 02" 


BL Lac 


0.028 


<0.08 


>950 


3 


PG 2209+184 


22^ 11™ 53.9= 


+18° 41' 50" 


RG 


0.070 


<0.13 


>400 


7 


RBS 1888 


4311141.6= 


-12° 31' 38" 


HBL 


0.226 


<0.009 


>170 


5 


HS 2250+1926 


221^ 53™ 07.4= 


+19° 42' 35" 


FSRQ 


0.284 


<0.009 


>590 


1 


2QZ J225453.2-272509 


221^ 54™ 53.2= 


-27° 25' 09" 


BL Lac 


0.333 


<0.016 


>170 


5 


PKS 2254+074 


221^ 57™ 17.3= 


+07° 43' 12" 


LBL 


0.193 


<0.05 


>900 


3 


NGC 7469 


23!^ 03™ 15.6= 


+08° 52' 26" 


Syl 


0.017 


<0.006 


>250 


5 


PKS 2316-423 


23!^ 19™ 05.8= 


-42° 06' 49" 


HBL 


0.055 


<0.014 


>190 


5 


lES 2321+419 


231^23™ 52.5= 


+42° 10' 55" 


HBL 


0.059 


<0.03 


>890 


3 


lES 2343-151 


23!^ 45™ 38.4= 


-14° 49' 29" 


IBL 


0.224 


<0.012 


>230 


1 



•^BL Lac, LBL, IBL, HBL - BL Lac object; FSRQ - Flat Spectrum Radio Quasar; FRl & 2 - Fanaroff- 
Riley 1 & 2 galaxy; Sy 1, 1.5 & 2 - Seyfert 1, 1.5 & 2; RG - Radio Galaxy; NLSyl - Narrow- 
Line Seyfert 1; UnC — Unclassified; see the SIMBAD l |http : //slmbad.u-stras bg. f r/simbad ) and NED 
(http : //ned www. ipac . caltech. edu[ | databases. 

''Tentative measurement, see comment in text. 
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Table 4. Target objects with associations in the OFGL source hst and 3EG catalog.. 



TeV name OFGL name 3EG name 



TeV detected: 



3C 66A 


OFGL 


J0222 

fj \J 


6+4302 


3EG 


10222+4253 


S5 0716+714 


OFGL 


T0722 


0-1-7120 


3EG 


70721+7120 


lES 1011+496 


OFGL 


J1015 


2+4927 






Markarian 421 


OFGL 


J1104 


5+3811 


3EG 


Jl 104+3809^ 


W Comae 


OFGL 


J1221 


7+2814 


3EG 


J1222+2841 


3C 279 


OFGL 


J1256 


1-0547 


3EG 


J1255-0549 


PKS 1424+240 


OFGL 


J1427 


1+2347 






PG 1553+113 


OFGL 


J1555 


8+1110 






Markarian 501 


OFGL 


J1653 


9+3946 






lES 1959+650 


OFGL 


J2000 


2+6506 






PKS 2005-489 


OFGL 


J2009 


4-4850 






PKS 2155-304 


OFGL 


J2158 


8-3014 


3EG 


J2158-3023^ 


BL Lacertae 


OFGL 


J2202 




3EG 


J2202+4217 


Centaurus A 


OFGL 


J1325 


4-4303 


3EG 


J1324-4314 


TeV non-detected: 










lES 0033+595 


OFGL 


J0036 


7+5951 






ON 325 


OFGL 


J1218 


0+3006 






3C 273 


OFGL 


J1229 


1+0202 


3EG 


J1229+0210 


NGC 1275 


OFGL 


J0320 


0+4131 







well established TeV source at the time of the 3EG 
catalog. 
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Table 5. Fermi-LAT detections (0.2 GeV-300 GeV). 







Parameters of fitted power-law spectrum 




Highest 


energy 


ProbabiUty of 


Name 


TS 


r iux (>zUUMeVj 




Photon Index 


Decorr. 


photons 


constant flux 






TP -L. A IT 1 A T-i 


r ± APstat ± AFsys 


energy 


1 St 


cth 



iU day 


zo uay 




[1] 


[iU cm s J 












[GeV] 




[oe vj 


ri 1 


rii 
lij 


TeV detected: 


























3C 66A 


2221 


96.7 ± 5.82 ± 3.39 


1 


.93 


± 0.04 


± 


0.04 


1.54 


111=^ 


54 


< 0.01 


< 0.01 


RGB J0710+591 


42 


0.087 ±0.049 ± 0.076 


1 


.21 


± 0.25 


± 


0.02 


15.29 


74 


4 


0.98 


0.94 


S5 0716+714 


1668 


79.9 ± 4.17 ± 2.84 


2, 


.16 


± 0.04 


± 


0.05 


0.82 


63 


9 


< 0.01 


< 0.01 


lES 0806+524 


102 


2.07 ± 0.38 ± 0.71 


2, 


.04 


± 0.14 


± 


0.03 


1.54 


30 


4 


0.05 


< 0.01 


lES 1011+496 


889 


32.0 ± 0.27 ± 0.29 


1 


.82 


± 0.05 


± 


0.03 


1.50 


168 


32 


0.54 


0.50 


Markarian 421 


3980 


94.3 ± 3.88 ± 2.60 


1 


.78 


± 0.03 


± 


0.04 


1.35 


801 


155 


0.06 


0.02 


Markarian 180 


50 


5.41 ± 1.69 ± 0.91 


1 


.91 


± 0.18 


± 


0.09 


1.95 


14 


2 


0.98 


0.54 


lES 1218+304 


147 


7.56 ± 2.16 ± 0.67 


1 


.63 


± 0.12 


± 


0.04 


5.17 


356 


31 


0.53 


0.06 


W Comae 


754 


41.7 ± 3.40 ± 2.46 


2, 


.02 


± 0.06 


± 


0.05 


1.13 


26 


18 


0.01 


< 0.01 


3C 279 


6865 


287 ± 7.13 ± 10.2 


2, 


.34 


± 0.03 


± 


0.04 


0.59 


28 


21 


< 0.01 


< 0.01 


PKS 1424+240 


800 


34.35 ± 2.60 ± 1.37 


1 


.85 


± 0.05 


± 


0.04 


1.50 


137 


30 


< 0.01 


0.16 


H 1426+428 


38 


1.56 ± 1.05 ± 0.29 


1 


.47 


± 0.30 


± 


0.11 


8.33 


19 


3 


0.83 


0.39 


PG 1553+113 


2009 


54.8 ± 3.63 ± 0.85 


1 


.69 


± 0.04 


± 


0.04 


2.32 


157 


76 


0.40 


0.54 


Markarian 501 


649 


22.4 ± 2.52 ± 0.13 


1 


.73 


± 0.06 


± 


0.04 


2.22 


127 


50 


0.57 


0.18 


lES 1959+650 


306 


25.1 ± 3.49 ± 2.83 


1 


.99 


± 0.09 


± 


0.07 


1.60 


75 


21 


0.91 


0.29 


PKS 2005-489 


246 


22.3 ± 3.09 ± 2.14 


1 


.91 


± 0.09 


± 


0.08 


1.01 


71 


8 


0.86 


0.97 


PKS 2155-304 


3354 


109 ± 4.45 ± 3.18 


1 


.87 


± 0.03 


± 


0.04 


1.13 


299 


46 


< 0.01 


< 0.01 


BL Lacertae 


310 


51.6 ± 5.81 ± 12.2 


2, 


.43 


± 0.10 


± 


0.08 


0.85 


70 


4 


0.61 


0.23 


lES 2344+514 


37 


3.67 ± 2.35 ± 1.62 


1 


.76 


± 0.27 


± 


0.23 


5.28 


53 


3 


0.76 


0.46 


M 87 


31 


7.56 ± 2.70 ± 2.24 


2, 


.30 


± 0.26 


± 


0.14 


1.11 


8 


1 


0.43 


0.57 


Centaurus A 


308 


70.8 ± 5.97 ± 5.80 


2, 


.90 


± 0.11 


± 


0.07 


0.47 


6 


4 


0.38 


0.97 


TeV non-detected: 


























lES 0033+595 


137 


20.3 ± 5.11 ± 1.74 


2, 


.00 


± 0.13 


± 


0.07 


2.58 


150 


16 


0.40 


0.01 


MG J0509+0541 


217 


19.7 ± 3.78 ± 0.70 


2, 


.01 


± 0.11 


± 


0.06 


1.95 


31 


12 


0.73 


0.23 




1 zlS 
14(5 


26.6 ± 3.50 ± 3.34 


2, 


.52 


± 0.13 


± 


0.10 




7 


2 


0.03 


0.11 


lES 0647+250 


95 


4.09 ± 1.39 ± 1.01 


1, 


.66 


± 0.15 


± 


0.09 


4.54 


247 


16 


0.30 


0.72 


PKS 0829+046 


187 


27.3 ± 3.37 ± 1.08 


2, 


.43 


± 0.11 


± 


0.04 


0.70 


4 


2 


0.38 


0.11 


lES 1028+511 


52 


3.88 ± 1.43 ± 0.57 


1, 


.72 


± 0.19 


± 


0.08 


3.07 


48 


2 


0.85 


0.31 


RGB J1117+202 


116 


7.12 ± 1.75 ± 0.36 


1, 


.79 


± 0.13 


± 


0.06 


2.37 


46 


6 


0.63 


0.97 


lES 1118+424 


33 


2.31 ± 1.27 ± 0.41 


1, 


.71 


± 0.26 


± 


0.08 


3.99 


27 


3 


0.64 


0.69 


ON 325 


761 


42.3 ± 3.68 ± 2.98 


1, 


.99 


± 0.06 


± 


0.06 


1.26 


45 


12 


< 0.01 


< 0.01 


3C 273 


3569 


224. ± 6.78 ± 8.49 


2, 


.79 


± 0.04 


± 


0.04 


0.45 


11 


5 


< 0.01 


< 0.01 


RX J1417+2543 


31 


2.56 ± 2.14 ± 0.65 


1, 


.68 


± 0.39 


± 


0.08 


6.02 


41 


1 


0.95 


0.32 


lES 1440+122 


33 


1.05 ± 0.06 ± 0.10 


1, 


.18 


± 0.27 


± 


0.03 


17.04 


19 


2 


0.68 


0.86 


RX J1725.0+1152 


152 


18.1 ± 3.93 ± 0.86 


2, 


.01 


± 0.13 


± 


0.05 


1.87 


39 


11 


< 0.01 


0.02 


I Zw 187 


31 


5.41 ± 2.23 ± 0.68 


1, 


.95 


± 0.23 


± 


0.03 


2.22 


77a 


3 


0.78 


0.41 


lES 1741+196 


46 


4.93 ± 2.17 ± 0.17 


1, 


.80 


± 0.22 


± 


0.03 


3.58 


37 


3 


0.81 


0.99 


lES 2321+419 


88 


6.76 ± 2.77 ± 0.79 


1, 


.78 


± 0.20 


± 


0.09 


3.90 


42 


14 


0.34 


0.13 


NGC 1275 


1351 


99.1 ± 5.13 ± 3.87 


2, 


.20 


± 0.04 


± 


0.06 


0.80 


18 


13 


< 0.01 


< 0.01 



''Photon-like events were selected for study using the so-called "diffu se" class cuts. Using a stricter set of cuts ( "extradiffuse" ) , 
developed to study the extragalactic diffuse radiation dAbdo et al.ll2009eh . the highest energy photons from 3C 66A and I Zw 187 were 
eliminated, giving Emax = 90 GeV and 4 GeV respectively for the two sources. 
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Table 6. Parameters of fitted power-law spectra in low-energy (0.2 GeV-1 GeV) and 
high-energy bands (1 GeV-300 GeV). Only Fermi sources detected with TS > 100 in each 

band are listed. 



Low-energy band (0.2 GeV-1 GeV) High-energy band (1 GeV-300 GeV) 

Name Flux Photon Index Decorr. Flux Photon Index Decorr. 





F± AFstat 
[10-9cm-2s-i] 


r ± AFstat 
[1] 


energy 
[GeV] 


F± AFstat 
[10-9cm-2s-i] 


F ± AFstat 
[1] 


energy 
[GeV] 


TeV detected: 














30 66A 


80.5 ± 7.2 


1.97 ± 0.16 


0.52 


17.3 ± 1.0 


1.98 ± 0.04 


2.47 


S5 0716+714 


65.8 ± 1.7 


2.20 ± 0.05 


0.37 


12.8 ± 2.2 


2. ,37 ± 0.09 


2.17 


lES 1011+496 


23.2 ± 3.1 


2.11 ± 0.25 


0.47 


4.6 ± 0.8 


1.96 ± 0.09 


3.08 


Markarian 421 


74.4 ± 4.2 


1.93 ± 0.11 


0.47 


9.6 ± 0.9 


1.78 ± 0.04 


3.71 


W Comae 


33.7 ± 3.9 


1.92 ± 0.21 


0.51 


5.5 ± 1.1 


2.16 ± 0.10 


2.52 


30 279=- 


142 ± 7.61 


2.49 ± 0.11 


0.40 


16.2 ± 1.5 


2.55 ± 0.12 


1.98 


30 279'' 


512.6 ± 18.3 


2.00 ± 0.08 


0.43 


67.9 ± 4.5 


2.49 ± 0.09 


1.98 


PKS 1424+240 


25.9 ± 3.0 


1.84 ± 0.22 


0.50 


8.41 ± 0.63 


1.82 ± 0.05 


1.77 


PG 1553+113 


34.8 ± 4.6 


1.52 ± 0.25 


0.58 


5.6 ± 0.6 


1.70 ± 0.05 


4.14 


PKS 2155-304 


78.6 ± 4.5 


1.72 ± 0.11 


0.49 


13.2 ± 1.2 


1.96 ± 0.04 


2.71 


TeV non-detected: 














ON 325 


32.3 ± 4.4 


1.98 ± 0.23 


0.48 


8.5 ± 1.7 


2.32 ± 0.11 


2.25 


30 273 


206.5 ± 6.8 


2.66 ± 0.07 


0.38 


62.5 ± 17.7 


3.42 ± 0.17 


1.47 


NGC 1275 


81.6 ± 5.4 


2.09 ± 0.13 


0.46 


14.5 ± 2.9 


2.40 ± 0.11 


2.25 



■^Pre-flaring period (MJD< 54780). 
''Peak flaring (54790 <MJD< 54830). 



Table 7. Fermi-LAT 95% flux upper-limits (0.2 GeV-300 GeV) assuming spectral indices 

of r = 1.5 and T = 2.O.. 



Flux limit, assuming 
Name T = 1.5 T = 2.0 

[lO-'^cm-^s-^] 



TeV detected: 



RGB J0152+017 


2.02 


5.01 


lES 0229+200 


1.94 


5.12 


lES 0347-121 


0.80 


1.81 


PKS 0548-322 


0.59 


3.14 


lES 1101-232 


0.83 


4.40 


H 2356-309 


0.28 


7.25 
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Table 8. Extrapolation of measured GeV spectrum into TeV regime. 



Extrapolation at 100 GeV Extrapolation over 0.2 TeV to 10 TeV band 
Name dF/dE (100 GeY) Photon index - Text Integral flux - </)ext 

[10-9cm-2s-lTeV-l] [1] [^Crab] 



TeV detected: 



3C 66A 


1.98 


± 


0.34 


4.26 


± 


0.02 


0.0750 ± 0.0142 


RGB J0710+591 


0.61 


± 


0.38 


1.74 


± 


0.11 


0.2110 + 0.2015 


85 0716+714 


0.34 


± 


0.12 


3.20 


± 


0.03 


0.0135 ± 0.0060 


lES 0806+524 


0.18 


± 


0.11 


2.65 


± 


0.10 


0.0184 ± 0.0143 


lES 1011+496 


1.05 


± 


0.33 


2.82 


± 


0.04 


0.0913 ± 0.0368 


Markarian 421 


5.86 




0.92 


1.90 


± 


0.02 


1.4351 ± 0.3265 


Markarian 180 


0.17 




0.12 


2.09 


± 


0.13 


0.0315 ± 0.0307 


lES 1218+304 


0.94 


± 


0.38 


2.46 


± 


0.07 


0.1354 ± 0.0739 


"W Comae 


0.46 


± 


0.18 


2.44 


± 


0.04 


0.0483 ± 0.0251 


3C 279 


0.35 


± 


0.09 


5.25 


± 


0.02 


0.0058 ± 0.0017 


PKS 1424+240 


1.47 


± 


0.30 


1.85 


± 


0.03 


0.4187 + 0.1236 


H 1426+428 


0.42 


± 


0.32 


1.85 


± 


0.18 


0.0885 ± 0.0108 


PG 1553+113 


3.32 


± 


0.56 


5.27^ 


± 0.02 


0.0423 ± 0.0074 










2.06'' 


± 0.03 


0.9960 ± 0.2327 


Markarian 501 


1.71 


± 


0.44 


1.86 


± 


0.05 


0.4414 ± 0.1572 


lES 1959+650 


0.53 


± 


0.20 


2.17 


± 


0.06 


0.0867 ± 0.0433 


PKS 2005-489 


0.68 


± 


0.26 


2.20 


± 


0.06 


0.1121 ± 0.0555 


PKS 2155-304 


3.11 




0.54 


2.37 


± 


0.02 


0.3942 ± 0.0882 


BL Laccrtac 


0.10 


± 


0.05 


2.72 


± 


0.08 


0.0081 ± 0.0048 


lES 2344+514 


0.23 


+ 


0.20 


1.94 


± 


0.16 


0.0542 ± 0.0690 


Centaurus A 


0.010 


± 


0.006 


2.90 


± 


0.11 


0.0006 ± 0.0004 


M 87 


0.027 


± 


0.033 


2.33 


± 


0.22 


0.0033 ± 0.0052 


TeV non-detected: 
















lES 0033+595 


0.39 


± 


0.19 


2.35 


± 


0.08 


0.0532 ± 0.0351 


MG J0509+0541 


0.36 


± 


0.17 


2.01 


± 


0.07 


0.0708 ± 0.0446 


PKS B0521-365 


0.03 


± 


0.02 


2.74 


± 


0.14 


0.0025 ± 0.0021 


lES 0647+250 


0.43 


± 


0.22 


2.61 


± 


0.08 


0.0543 ± 0.0363 


PKS 0829+046 


0.05 


± 


0.03 


3.24 


± 


0.09 


0.0029 ± 0.0020 


lES 1028+511 


0.28 




0.20 


3.56 


± 


0.10 


0.0180 + 0.0162 


RGB ,11117+202 


0.40 




0.22 


2.40 


± 


0.09 


0.0572 ± 0.0398 


lES 1118+424 


0.19 


± 


0.17 


2.25 


± 


0.15 


0.0319 ± 0.0391 


ON 325 


0.28 


± 


0.12 


3.12 


± 


0.04 


0.0151 ± 0.0079 


3C 273 


0.05 




0.01 


3.52 


± 


0.03 


0.0001 ± 0.0001 


RX J1417+2543 


0.24 


± 


0.28 


2.79 


± 


0.18 


0.0262 ± 0.0413 


lES 1440+122 


0.80 


± 


0.48 


1.68 


± 


0.11 


0.2310 ± 0.2178 


RXJ 1725.0+1152 


0.33 




0.18 


2.08 


± 


0.10 


0.0601 ± 0.0452 


I Zw 187 


0.14 


± 


0.13 


2.17 


± 


0.16 


0.0235 ± 0.0287 


lES 1741+196 


0.27 


± 


0.20 


2.14 


± 


0.13 


0.0484 + 0.0511 


lES 2321+419 


0.41 


± 


0.29 


2.01 


± 


0.13 


0.0868 ± 0.0865 


NGC 1275 


0.34 




0.14 


2.26 


± 


0.04 


0.0355 ± 0.0200 



''Extrapolated assuming z = 0.09. 
''Extrapolated assuming z = 0.78. 
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3C66A/B 1ES 1218+304 




1ES 1028+511 1ES 1118+424 




Fig. 1. — Selected TS maps, covering a 1° x 1° region around the source of interest. In each 
case, the location of the AGN is indicated with a "x". In the case of the 3C 66A/B field 
(top left), the location of 3C 66A is indicated with a "x" and that of 3C 66B with a "+". 
TS maps for all the targets in the study are available in the online materials. 
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10" 10'° 10" 10'^ ,rr,„10'^ 10° 10'° 10" 10"^ ,^r,o10'^ 

E [eV] E [eVl 



Fig. 2. — Spectra for the 21 GeV-TeV detected objects. The GeV spectrum derived from 
Fermi-LAT observations is indicated as a "butterfly" contour (sohd hne). For brighter 
sources (those in Table [6]) the contours correspond to the high-energy band [E > IGeV), 
with the fits over the full energy range shown as gray bands. For the weaker sources, only the 
fits over the full range (given in Table [5]) are shown. TeV spectral measurements published 
by H.E.S.S. (circles), VERITAS/Whipple (squares) and MAGIC (triangles) are also shown. 
An extrapolation of the Fermi spectrum to the TeV regime is shown (dashed line), assuming 
absorption with the EBL as described in the text. In the panel for the 3C 66A/B region, 
the extrapolation is shown for z = 0.444 (3C 66A - dashed line) and z = 0.021 (3C 66B - 
dash-dotted line). In the case of PG 1553+113 extrapolations with z = 0.78 (dashed line) 
and z = 0.09 (dash-dotted line) are shown. 
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Fig. 2. — Continued 
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10° 10" 10" 10 



E[eV] 



10" 10'" 10" 10'^ ^, ,J0 

E[eV] 



PKS 2155-304 



10' 10'° 10" 10'= ^, ,„10" 

EEeV] 



BLLac 




10' 10'° 10" io'= ,_io" 

E[eV] 




M87 




10° 10'° 10" io'= ,_io" 

E[eV] 



10° 10" 10" 10 



E[eV] 



Fig. 2. — Continued 
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Fig. 3. — Fermi upper-limits with a spectral index frozen at 1.5 and 2 for the TeV detected 
sources. See caption of Figure [2] for details. 
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Fig. 4. — 28-day light curves for Fermi-detected sources, centered on tlie new moon. For 
each flux point, the vertical bar shows the statistical error (only) and the horizontal bar 
indicates the duration of integration. For each source, the mean flux over the full duration 
of the study is shown as a dashed line, and the systematic uncertainty in the flux points, 
which is estimated to be 3% of the flux, is shown as a gray band. 
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Fig. 4. — Continued 
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Fig. 4. — Continued 



Date - MJD-54000 [days] 
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Fig. 5. — 10-day light curves for selected Fermz-detected sources. See caption of Figure 
for explanation of what is indicated on each panel. 
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Fig. 6. — Difference, AF, between the measured TeV and Fermi photon indices as a function 
of the redshift. Empty circles denote the BL Lacs, the filled circle denotes Cen A and the 
filled square M 87. 



